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ABSTRACT 


It  is  the  purpose  of  this  research  study  to  investigate  theorectical ly  several 
microwave  circuits  which  are  suitable  to  integrate  with  a  beam-plasma  millimeter 
wave  tube,  and  the  problems  associated  with  them. 

Three  microwave  structures  have  been  studied:  (1)  The  hollow  dielectric  tube 
waveguides,  (2)  The  cylindrical  metal  waveguide  partially  filled  with  a  coaxial 
dielectric  tube,  and  (3)  The  hollow  dielectric  tube  waveguide  supported  by  a 
conducting  cylinder.  They  all  employ  dielectric  tube  as  their  basic  element,  and 
they  all  operate  in  the  circularly  symmetric  transverse  electric  modes.  The 
field  problems  have  been  solved  and  the  resulting  transcendental  mode  equations 
are  solved  graphically  in  conjunction  with  their  respective  despersion  relations. 
It  is  found  that  a  critical  frequency  exists  for  all  three  structures,  and  below 
which  no  propagation  waves  can  take  place.  As  a  circuit  element  or  cavity  both 
types  I  and  II  structures  show  to  have  the  very  desirable  property  of  confining, 
and  re-directing  the  emitted  electromagnetic  radiation  from  the  plasma  tube. 

More  than  that  both  these  structures  are  capable  of  single  mode  operation  at 
very  short  wavelengths. 

It  is  a  surprise  to  find  that  when  type  II  structure  is  operated  at  the  very 
high  frequency  mode,  the  structure  behaves  as  if  it  is  a  dielectric  tube  wave¬ 
guide,  limited  on  the  outside  by  a  conducting  surface  of  the  cylindrical  metal 
wall,  and  on  the  inside  by  a  medium  of  small  dielectric  constant.  The  fields 
and  waves  are  confined  to  the  space  occupied  by  the  high  dielectric  constant 
dielectric  tube. 

Design  information  regarding  to  the  choice  of  the  waveguide  dimensions  and 
the  dielectric  constant  for  the  dielectric  tube  in  relation  with  the  guide 
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I.  Introduction 


One  of  the  more  pressing  problems  in  present  day  microwave  technology  is 
that  of  producing  a  tunable  source  of  coherent  radiation  in  the  millimeter  and 
submillimeter  wavelengths  (from  1  to  0.1  mm).  The  potential  use  of  millimeter 
waves  in  engineering  applications  are  almost  too  numerous  to  mention.  For 
example,  in  radar  system  design  millimeter  waves  could  give  us  high  resolution 
radar  for  guided  missiles,  precision  gun  and  rocket  fire  control,  accurate 
stationkeeping  in  fleet  formations...,  to  name  just  a  few.  For  communication 
engineers  we  like  to  think  of  the  progress  in  millimeter  wave  field  as  extend¬ 
ing  the  frontiers  of  communication  electronics  beyond  the  crowded  microwave 
spectrum,  and  into  a  new  frequency  region  where  the  spectrum  density  is  low 
and  signal  channels  less  crowded.  The  very  high  frequencies  and  the  quasi- 
optical  properties  which  characterize  millimeter  and  submillimeter  waves  would 
enhance  many  secret  conmuni cation  systems,  both  within  the  atmosphere  and  above 
it.  Very  long-range  pencil-beam  communications  may  be  easily  realized  and  yet 
require  only  a  small  antenna.  These  and  many  other  possible  applications  to 
enginering  system  designs  and  military  systems,  particularly  communications 
and  navigations,  have  long  been  recognized.  Thus,  from  the  technological 
point  of  view,  the  availability  of  a  coherent  tunable  millimeter  wave  power 
source  would  make  possible  a  significant  advance  in  all  these  and  many  other 
special  purpose  conmuni cation  systems. 

The  biggest  road  block  In  the  full  realization  on  these  advantages  is 
the  engineering  one  of  constructing  an  efficient  coherent  source  of  millimeter 
wave  tubes  which  Is  capable  of  providing  adequate  power  at  millimeter  wave¬ 
lengths  and  below.  The  various  problems  associated  with  fragile  RF  structures 
and  high  density  electron  beams  have,  up  to  the  present  time,  blocked  the 


extension  of  conventional  microwave-tube  techniques  down  to  one-millimeter  or 
shorter  wavelengths.  At  these  very  short  wavelengths,  conventional  microwave 
tubes  would  be  so  tiny  that  watch  making  techniques  are  required.  Furthermore, 
in  conventional  microwave  and  millimeter  wave  devices,  beam  interaction  and  RF 
power  dissipation  on  delicate  interaction  structures  limit  power  handling 
capacity. 

The  millimeter  wave  tube  described  in  the  next  section  is  a  beam-plasma 
device  based  on  the  electron  beam  plasma  interaction  in  warm  plasma.  The 
radiation  frequency  of  the  device  is  not  determined  by  a  resonant  circuit 
(e.g.,  the  cavity  or  the  slow  wave  structure  used  in  klystron  or  traveling  wave 
tube),  but  by  the  plasma  frequency  of  oscillation.  The  idea  of  utilizing  a 
plasma  oscillation  frequency  as  a  parameter  to  determine  the  radiation  fre¬ 
quency  of  the  device  is  an  attractive  one.  The  frequency  of  electron  oscill¬ 
ation  in  a  plasma  is  the  nature  frequency  associated  with  the  number  density 
of  the  electrons  in  plasma  which  can  be  controlled  by  gas  discharge;  thus  no 
external  radio  frequency  circuit  is  required  for  the  determination  of  its 
radiation  frequency. 

The  beam-plasma  device  offers  other  attractions  too.  For  example,  in  a 
beam-plasma  device,  high  power  can  be  handled  at  high  frequencies  since  the 
delicate  mechanical  structures  such  as  the  slow  wave  structures  in  the  travel¬ 
ing-wave-tube  or  the  cavity  resonators  in  the  klystron  are  replaced  by  the 
gaseous  plasma  volume  which  can  be  relatively  huge  in  physical  size.  Besides, 
a  large  size  electron  beam,  which  has  been  on  of  the  bottlenecks  in  extending 
the  conventional  microwave  tube  into  millimeter  wave  region,  can  be  employed 
prolifically  to  enhance  the  device  output  power. 

Although  the  power  output  of  the  beam-plasma  millimeter  wave  device  depends 
critically  upon  the  Injected  beam  current  and  beam  power,  only  weakly  relati- 


vistic  beam  may  be  needed  for  its  satisfactory  working.  When  this  is  compared 
with  the  very  large  relativistic  beam  used  in  NRL  megawatt  millimeter  wave  tube^2! 
the  beam-plasma  millimeter  wave  device  shows  greater  potentiality  as  a  high  power 
millimeter-wave  device.  Theoretically,  no  dc  magnetic  field  is  required  for  its 
operation;  however,  small  externally  applied  dc  magnetic  field  may  be  required 
for  the  confinement  of  the  electron  beam  motion  in  the  beam-plasma  millimeter 
tube,  but  this  is  to  be  contrasted  with  the  super-conducting  magnet  needed  for 
their  operation  in  a  new  class  of  millimeter  wave  tubes ^  employing  cyclotron 
fast  wave  interaction  such  as  the  gyratrons. 


II.  A  BRIEF  REVIEW  OF  THE  PRASAD-LEIBYS 1 2 3 4 5 6  BEAM-PLASMA  DEVICES 


The  physical  principle  and  gain  mechanism  of  the  becm-plasma  device  to  be 
described  in  the  following  are  based  on  Prasad-Leibys '  two  patents ^  and 
Prasad's  theoretical  paper^.  Basically  Prasad-Leibys'  plasma  tube  is  a 
power  oscillator  consisting  of,  in  its  most  elementary  form,  a  beam-plasma 
interaction  chamber  and  a  source  of  directed  un-modulated  high  velocity  but 
non-relativistic  electrom  as  sketched  in  fig.  1. 


1  -  Electron  gun 

2  -  Electron  beam 

3  -  Metallic  foil 

4  -  Plasma  chamber  (method  of  plasma 

generation  now  shown) 

5  -  Collector 

6  -  Cylindrical  glass  tube 

Figure  1.  Sketch  of  the  Beam-Plasma  Device. 


The  plasma,  which  may  be  either  discharge  generated  or  electron  beam  gen¬ 
erated,  is  unmagnetized,  warm  uniform,  and  collisionless.  The  electron  beam 
which  may  be  produced  by  any  one  of  the  electron  gun  systems  employed  in  kly¬ 
strons  and  traveling-wave  tubes  is  of  the  low  density  type,  having  a  density 
(nb)  about  2  to  5%  of  the  plasma  electron  density  (ng). 


J 


An  electron  beam  traversing  through  a  warm  plasma  excites  electron 
oscillations  in  the  plasma.  The  waves  associated  with  such  oscillations  are 
the  Langmuir  electron  waves  which  are  non-radiative  longitudinal  electrostatic 
waves.  They  cannot  escape  from  the  plasma  volume.  The  plasma  electrons  are 
accelerated  by  the  waves  propagating  through  it.  The  accelerated  electrons 
in  turn  emit  radiation  in  both  the  transverse  and  longitudinal  directions. 

If  the  plasma  was  completely  uniform,  the  radiation  emitted  would  simply 
reproduce  the  waves  themselves.  However,  in  a  warm  plasma,  the  presence  of 
both  ion  and  electron  density  fluctuations  result  in  a  scattering  and  coupling 
of  various  types  of  waves. 

Since  electrons  tend  to  follow  the  ions  so  as  to  maintain  charge  neutrality, 
there  will  be  electron  density  fluctuations  associated  with  ion  density  fluctu¬ 
ations.  There  are  also  density  fluctuations  associated  with  Langmuir  electron 
oscillations,  as  well  as  density  fluctuations  due  to  the  random  motion  of  the 
electrons.  The  scattering  of  these  waves  on  the  inhomogeneities  of  the  plasma 
ion  and  electron  distribution  results  in  the  generation  of  transverse  electro¬ 
magnetic  radiation  at  twice  the  plasma  frequency. 

The  possibility  of  transforming  the  energy  of  longitudinal  plasma  waves 
into  electromagnetic  radiation  in  a  homogeneous  and  fully  ionized  hot  plasma 
was  first  pointed  out  by  Ginburg  and  Zhelezhiakov^ .  In  essence,  the  gain 
mechanism  and  the  generation  of  transverse  electromagnetic  radiation  of  the 
Prasad-Leibyo'  millimeter  wave  tube  involves  a  three  step  beam- plasma  inter¬ 
action  process, 

(a)  Excitation  of  longitudinal  electrostatic  (or  pump)  waves  in  the  plasma 
by  the  injection  of  fast  electron  beam. 
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(b)  Generation  of  the  background  induced  plasma  waves  via  the  "Rayleigh 
Scattering"  of  the  pump  waves  on  the  inhomogeneities  of  the  plasma 
ion  distribution. 

(c)  Generation  of  the  transverse  electromagnetic  radiation  via  the  "Raman 
Scattering"  of  the  induced  background  plasma  waves  on  the  inhomogen¬ 
eities  in  the  plasma  electron  distribution. 

When  a  warm,  uniform,  collisionless  plasma  is  traversed  by  a  beam  of  fast, 
but  non-relativistic  charged  particles,  spatially  growing  waves  will  result 
from  the  excitation  of  oscillation  of  plasma  electron  by  the  beam  through  a  non¬ 
linear  interaction  mechanism,  and  the  interaction  of  oscillating  plasma  electron 
back  to  the  beam  (ion-acoustive  wave  will  be  excited  too).  The  wave  amplitude 
increases  exponentially  with  distance  along  the  beam,  and  the  rate  of  growth 
for  a  constant  excitation  frequency  rises  sharply  near  the  plasma  frequency, 
cOpg-  In  the  meantime,  the  ion-acoustic  waves  act  as  a  trigger  for  the  in¬ 
stability  and  lock  the  initial  phase  of  the  unstable  waves.  Moreover,  if  the 
instability  is  driven  by  a  cold,  weak  electron  beam  in  a  plasma,  it  will  be 
stabilized  by  the  trapping  of  plasma  particles  mechanism.  The  spectrum  of 
these  unstable  waves  is  so  narrow  that  only  a  single  frequency  grows;  in  the 
vicinity  of  the  point  where  the  wave  is  excited  in  the  plasma,  the  particles 
are  trapped  in  the  wave  potential  wells.  The  trapping  of  the  plasma  particles 
limits  the  instability,  and  the  trapped  particles  cause  oscillation  in  the 
amplitude  i  '  the  excited  wave  and  the  growth  of  the  new  waves,  the  so-called 
sidebands.  3 hi s  trapping  effect  has  been  analyzed^  and  experimentally 
observed^.  Such  waves  are  narrow  pockets  of  the  quasi-stable  Langmuir  waves^10^ 
which  propagate  in  the  direction  of  the  electron  beam.  The  wave  number  and 


frequency  distribution  are  peaked^^  at  (kQ,  wQ)  with  w0^wpe«  and  k0  =  wQ/vb 

(12) 

The  wave  number  kQ  and  frequency  satisfy  the  dispersion  relation' 

=  CO*  +  3  kn  vtP  .  (1) 

o  pe  o  te 

where 

<Jpe  ’  <4  *  ne  e2  '  "e  >  * 

_3 

is  the  electron  plasma  frequency,  with  ng  the  average  electron  density  in  cm  , 
e  charge  density  and  mg  the  mass  of  electron,  and 

»te  -  <k  Te  /  »e>  * 

is  the  thermal  velocity  of  the  electron,  with  Tg  the  electron  temperature,  k. 
Boltzmann  Constant. 

Such  waves  can  propagate  in  the  plasma  with  practically  no  (Landau) 
dumping  if  ajQ  is  close  to  oDpe  (i.e.,  when  kQ  is  very  small  compared  with 
the  inverse  Debye  Length  (  k  «  /\ ^*  =  /v.  ).  Furthermore,  the  condition 

for  propagation  is  that  the  phase  velocity  of  the  waves  (vph)  must  be  greater 
than  the  thermal  velocity  of  the  plasma  electron  (vte);  i.e.,  vpb  (crvb)  >  vtg, 
where  vb  is  the  beam  electron  drift  velocity,  and  for  undamped  wave  [.ropagation 
the  following  condition  must  satisfy  also: 

vph  <  *  vb  >  >  *te  .  (2) 

In  the  next  step,  as  the  incident  or  the  pump  waves  propagate  through  the 
plasma  medium,  they  are  scattered  by  the  ion  density  fluctuations  in  the  plasma 
The  scattered  wave  consists  of  two  component  waves:  (a)  the  strong,  but,  non- 
radiative  longitudinal  electrostatic  waves,  and  (b)  the  much  weaker  transverse 
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electromagnetic  wave  which  radiates  out  of  the  plasma  volume,  at  the  frequency 


The  frequency  of  the  scattered  wave  remains  nearly  the  same  as  that  of 
the  incident  pump  wave.  This  is  also  in  accordance  with  the  dispersion  rela- 
t1on<’2>; 


where  and  vtg  have  been  defined  in  Eq.  (1),  dJj  is  the  Rayleigh  scattered 
longitudinal  wave  and  k1  its  wave  number. 

This  is  known  as  the  "Rayleigh  Scattering"  process  in  plasma.  Now,  since 
the  Rayleigh  scattered  longitudinal  component  wave  is  incapable  of  escaping 
from  the  plasma  volume,  it  will  raise  the  amplitude  of  space  charge  density 
fluctuation  well  above  the  thermal  level  in  the  bulk  of  the  plasma,  and  hence 
contributes  to  a  growing  coherent  background  induced  wavefield,  EL  (cO^k^) 
with  o)|  a)o  and  k-j  ~  kQ  in  the  pumped  but  stable  plasma.  The  angular  dis¬ 
tribution  of  the  scattered  wave  is  that  of  an  electric  dipole^  and  the 
spatial  distribution  of  the  Rayleigh  scattered  wave  power  is  shown  in  Fig.  2. 
This  is  a  polar  diagram.  The  k^  vector  is  in  the  direction  of  the  injected 
electron  beam.  We  notice  that  the  very  strong  forward-lobe  electric  field 
varies  greatly  as  the  ratio  k^/kQ  is  varied,  whereas,  the  much  weaker  back¬ 
ward-lobe  field  remain  practically  unaltered  for  a  similar  change  in  the  ratio 


of  kl/k0. 
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Ep  -  The  Forward-Lobe  Electric  Field 

Eg  -  The  Backward-Lobe  Electric  Field 

Figure  2.  Asymmetrical  Dipole  Spatial  Distribution  of  the  Rayleigh  Scattered 
Wave  Power. 

When  the  Rayleigh  Scattered  longitudinal  component  wave  which  is  the  non- 
radiative  induced  background  waves  (&},  ,k^),  propagate  through  the  plasma 
volume,  a  small  percentage  of  these  waves  are  backscattered  on  the  inhomogen¬ 
eities  in  the  charge  distribution  of  the  plasma  resulting  in  a  Raman  scattered 
transverse  electromagnetic  waves  which  radiate  out  of  the  plasma  volume  at 
twice  the  electron  plasma  frequency  ('Wo?  ).  Since  the  inhomogeneities,  or 
the  fluctuations  of  plasma  electron  density,  are  associated  with  the  Langmuir 
electron  wave,  which  in  turn,  producing  the  pump-waves.  It  appears  as  though 
this  scattering  and  generation  of  second  harmonic  waves  result  when  the  induced 
background  plasma  waves  scatter  off  their  cousins  -  the  waves  which  have  been 
generated  a  moment  before  by  some  preceding  particles  much  like  the  wakes  of 
a  long  line  of  speed  boats.  There  will  be  a  continual  overtaking  and  scatter¬ 
ing  of  one  wave  by  another  because  the  spectrum  of  the  induced  waves  is  dis¬ 
tributed  both  in  angle  and  in  frequency;  each  component  has  its  own  group 
velocity. 


Under  certain  conditions,  when  the  energy  of  the  incident  or  pump  wave 


exceeds  a  certain  threshold  condition,  the  waves  become  phase  locked  with 


each  other,  enhanced  Raman  emission  occurs,  and  the  resulting  radiation  be- 

(6) 


comes  coherent.  The  threshold  condition  is  given  by 

2 

o  o  m  n  v.  .  v.  6 

U*  >  6xl03  (  — )  (  Z"  )  (  #  ) 


(4) 


'ph 


where  UQ  is  the  amplitude  of  the  pump  wave  electric  field, 

v^  and  vte  are  the  thermal  ion  and  electron  velocity  respectively, 

nQ  is  the  average  background  plasma  electron  density, 

Vp^  and  fjpe  are  phase  velocity  and  electron  plasma  frequency  given 
before. 


This  coherence  of  the  electromagnetic  radiation  at  discrete  frequency  bands  by 
an  electron  beam  interacting  with  a  warm  plasma  have  been  experimentally  ob- 
served"').'13'. 

Since  the  Rayleigh  scattered  longitudinal  component  wave  has  a  superthermal 
phase  velocity,  [which  follows  from  Eq.  (3)]  which  is  comparable  to  that  of  the 
incident  pump  wave,  and  like  the  pump  waves  the  wave  is  undamped  (Landau),  and 
by  Gauss'  law  they  will  appear  as  the  charge  density  fluctuations,  V-  E ^ ( r cO ^ ) 
in  the  plasma  which  can  readily  interact  with  the  incident  pump  wave  field, 

EQ(r,  coQ);  also  interchangably,  the  field,  E^(r,G>j),  will  interact  with  the 
fluctuations  associated  with  the  pump  wave  field  through  the  divergence  relation, 
.  E^r,  6?o).  These  interaction  relations  appear  in  the  form 

S(r,t)  =  (e/m)  (E.  y  .  ^  ^  y  .  E. )  .  (5) 

which  is  the  source  function  Prasad  used  to  evaluate  the  Ra  in  scattered 
electric  field,  Ej  (r,  Moreover,  because  electron  density  fluctuates 

at  the  electron  plasma  frequency,  <^e,  which  is  comparable  to  the  frequency  of 


the  incident  pump  wave,  and  by  momentum  and  energy  conservation  principle,  the 

scattering  process  gives  rise  to  the  emission  of  enhanced  transverse  Raman 

(141  (151 

radiation  '  '  in  a  three-wave  interaction  process'  ,  that  is 


cOl 

- 

or  a)2 

"  *>1 

+ 

• 

-1 

-2 

‘  V 

or  k_2 

=  -1 

+  ^ 

.  (6) 

where  ^  and  ^  are  the  angular  frequencies  of  the  incident  pump,  Rayleigh 
scattered  longitudinal  electrostatic,  and  the  Raman  scattered  transverse  electro¬ 
magnetic  waves  respectively,  J^,  k-j  and  are  their  respective  wave  vectors. 

We  observe  that,  since  and  k^  are  much  greater  in  magnitude  than  k^ 
(except  for  relativistic  energy  beam),  transverse  waves  can  only  be  produced 
by  a  "head-on"  scattering  collisions  of  the  longitudinal  electrostatic  waves. 

This  is  in  accordance  with  the  three-wave,  or  coherent  interaction  process 
described  above.  These  are  radiative  transverse  electromagnetic  waves  which 
are  capable  of  escaping  from  the  plasma  volume.  Since  £•  cJQ  (£  a^e) 
it  follows  from  Eq.  (6)  that  a>2~  ■  CJpe.  The  radiation  frequency  is,  there¬ 
fore,  twice  the  electron  plasma  frequency.  This  process  of  non-linear  coupling 
between  two  electrostatic  waves  and  one  electromagnetic  wave  is  in  accordance 
with  the  mechanism  of  second  harmonic  generation  proposed  by  Sturrock  et  al^16^. 

In  the  far  field  approximation  the  electric  component  wave  of  the  Raman 
scattered  transverse  electromagnetic  wave  has  been  calculated^  and  it  is 


given  by: 

IT(r,  2) 

where 
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k2  *  ii  k2 


pe 


(7) 


and  n  =  r/r 


c 

O' 


is  the  unit  vector  specifying  the  arbitrary  direction  of 
observation  of  the  transverse  Raman  radiation  field, 
is  the  velocity  of  light. 

is  a  small  positive  real  number  signifying  an  outgoing  wave 
decreasing  (slowly)  in  amplitude  after  scattering  at  a  source 
situated  anywhere  inside  the  scattering  volume. 


The  angular  distribution  pattern  of  the  Raman  scattered  waves  is  that  of 
an  electric  quadrupole,  characteristic  of  a  spontaneous  Raman  scattering  pro¬ 
cess  in  plasma  ^)»  0$)  m  illustrated  in  Fig.  3. 


Figure  3.  The  directional  pattern  $(a,0)  for  the  emission  of  the  enhanced 
second-harmonic  radiatio.i  in  the  warm  uniform  plasma,  for  a*0.3, 
a=0.5,  and  a=0.6.  (By  Courtesy  of  B.  Prasad). 
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Since  the  Raman  scattered  wave  is  transverse  electromagnetic  waves,  the  man- 
netic  field  can  be  computed  readily  to  yield 

Bt  (r,  )  =  (  n  x  Et  )  |~1-  (cOpe/cA?2)2  "I  1  .  (8) 


The  time-average  energy  flux  of  Raman  scattered  transverse  electromagnetic 
radiation  can  be  computed,  and  the  volume  emissivity,  that  is,  the  total  inten¬ 
sity  per  unit  volume  determined,  the  result  is  given  by  ^ 


<2 


„  _ F{a}_ 

45  (2H ) 
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.  <0 o  vte 
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m  n. 


'te 


....  (9) 


where  the  factor  F  (a)  depends  upon  the  fluctuation  density  function  which 
characterized  the  plasma,  and  for  non-thermal  conditions  F  can  be  as  large 
as  10  to  100. 

The  maximum  interaction  power  density  for  the  Raman  scattered  radiation 

(41 

at  2a)  has  been  calculated  '  '  in  terms  of  physical  size  and  the  beam-plasma 
pe 

interaction  length  L,  and  the  injection  beam  current  density,  Jb,  as  well  as 
the  beam  power  P^.  The  result  obtained  is 

Pmax(2c0pe)  =  3.6819xl0-20  L  F .  (10) 

where  d  is  the  ratio  of  electrostatic  wave  phase  velocity,  vpb,  to  the  mean 

2 

thermal  velocity  of  the  plasma  electrons,  vtg,  L  is  in  cms,  in  amperes/cm  , 
and  Pb  in  watts/cc. 

15) 

In  their  second  patent'  ,  Prasad  and  Lei  by  used  counter-streaming  beams 
instead  of  a  single  electron  beam.  Based  on  their  computation,  resonant  Raman 
scattering  of  the  longitudinal  induced  waves  on  inhomogeneities  of  the  plasma 
electron  distribution  in  a  counter-streaming  electron-beam  plasma  system  is 
capable  of  generating  10  to  100  times  as  much  coherent  electromagnetic  radia¬ 
tion  (at~2  &)pe)  as  would  be  produced  by  a  single  beam-plasma  system.  In  a 


-  14  - 


single  beam-plasma  device,  generation  of  Raman  scattered  transverse  electro¬ 
magnetic  radiation  is  due  to  the  "head-on"  collisional  interaction  of  the 
Rayleigh-scattered  induced  background  waves  on  the  fluctuations  of  plasma  charge 
density,  hence,  only  a  small  portion  of  the  relatively  weak  Rayleigh-scattered 
backward-lobe  field,  Eg,  contributes  to  the  Raman  process.  The  much  stronger 
forward  scattered  lobe,  £p,  makes  little  or  no  contribution  to  this  process. 

It  is  clear  that  a  greatly  enhanced  coherent  Raman  transverse  radiation  output 
should  result  if  oppositely  directed  pump-waves  (such  as  would  be  excited  by 
counter- st reaming  electron  beam)  were  present  in  the  plasma  so  that  both  lobes 
of  both  sets  of  longitudinal  (Rayleigh-scattered)  induced  plasma  waves  would 
participate  in  the  Raman  scattering  process.  Schematically,  the  counter¬ 
streaming  beam-plasma  system  is  illustrated  in  fig.  4. 


■-  . .  .. 


....  .....  . 


1  -  Electron  gun 

2  -  Cylindrical  glass  tube 

3  -  Metallic  foil 

A  -  Right-directed  electron  beam 
B  -  Left-directed  electron  beam 

C  -  Beam-plasma  interaction  chamber  (method  of  plasma  production 
not  shown. 

Figure  4.  Schematic  Diagram  of  a  Counter-Streaming  Beam-Plasma  Device  Model 

The  physical  principle  of  operation  of  ciis  device  is  described  briefly  in 
the  following: 

If  the  source  function  of  the  single-beam  plasma  device  is  written  using 
double  super-  and  sub-  scripts  notation,  one  would  obtain  from  eq.  (5) 


where  the  super-script  "R"  refers  to  all  field  quantities  that  are  associated 
with  the  right-directed  electron  beam.  The  first  sub-script  in  the  two  electric 
field  vectors,  E^  and  E^  refer  to  the  Rayleigh-scattered  longitudinal  induced 
plasma  wave  field  and  the  pump  wave  field  respectively;  whereas  the  second  sub- 
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script  "B"  in  the  expression  of  denotes  the  field  that  belongs  to  the  back¬ 
ward  lobe  of  the  induced  plasma  wave  field.  The  double  super-and  sub-  scripts 
that  appear  in  the  source  function  S  are  used  to  identify  the  component  wave 

fields  that  are  directly  involved  in  the  "head-on"  collisional  interaction 

R  R 

producing  this  source  field.  Thus,  the  symbol  S  *  implies  that  this  source 

B.O 

function  has  its  origin  from  the  collisional  interaction  between  the  backward- 
lobe  electric  field  of  the  right-directed  electron  beam-generated  Rayleigh- 
scattered  longitudinal  induced  plasma  wave  and  the  pump  wave  field  excited 
by  the  right-directed  electron  beam. 

Using  this  symbolism,  we  can  write  down  the  source  function  that  has  its 
origin  from  the  "head-on"  collisional  interaction  between  the  backward-lobe 
electric  field  of  the  induced  plasma  wave  and  the  fluctuation  of  the  electron 
density  that  are  associated  with  a  left-directed  single-beam  as  follows: 


Accordingly,  the  angular  distribution  of  the  Rayleigh-scattered  induced 
background  wave  field  of  the  left -directed  single-beam  device  will  have  an 
asymmetric  dipole  distribution  like  that  of  fig.  2  except  that,  in  the  present 
case,  both  the  forward-  and  backward-lobes  have  rotated  by  180-degrees.  This 
is  illustrated  in  fig.  5. 


£  -  the  forward-lobe  electric  field 
F 

i  -  the  backward-lobe  electric  field 
B 

Figure  5.  Asymmetric  dipole  spatial  distribution  of  the  Rayleigh-scattered 
Induced  plasma  wave  field  of  a  left-directed  single  beam-plasma 
device. 
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The  source  functions  described  above  have  their  origins  from  the  inter¬ 
action  mechanism  arising  from  the  "head-on"  collision  between  the  pump  wave 
field  excited  by  one  electron  beam,  and  the  backward-lobe  electric  f!eld  „f 
the  induced  plasma  waves  which  have  been  generated  a  moment  before  by  some 
preceding  particles  of  the  same  electron  beam  -  much  like  the  wakes  of  a  long 
line  of  speed  boats.  In  a  counter-streaming  beam-plasma  device,  in  addition 
to  the  interaction  mechanisms  just  mentioned,  there  are  two  more  sets  of 
source  functions  arising  from  the  mutual  "head-on"  collisional  interaction 
between  the  pump  wave  field  excited  by  one  beam,  and  the  forward-lobe  electric 
field  of  the  induced  plasma  wave  field  generated  by  an  oppositely  directed 
beam.  They  are: 

(i)  Source  function  originated  from  the  "head-on"  collision  between  the 

p 

right-directed  electron  beam-excited  pump  wave  field  E  (w  ,k  )  and 

0 

the  forward-lobe  electric  field  of  the  Rayleigh-scattered  longitudinal 
induced  background  plasma  wave  field  F  (w^k^),  generated  by  the 
left-directed  electron  beam, 


that  is 
.L,R 


=  (e/m) 
F  ,0 


[  V 


er  +  er 


ii.F 


(13) 


(ii)  Source  function  originated  from  the  "head-on"  collision  between  the 

left-directed  electron  beam-excited  pump  wave  field,  E*"  (w  ,-k  ), 

0  0 

and  the  forward-lobe  electric  field  of  the  Rayleigh-scattered  longi- 

D 

tudinal  induced  plasma  wave,  E^  p  (w^.jcj),  generated  by  the  right- 
directed  electron  beam. 


that  is 
SR,L 
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Eqs.  (11),  (12),  (13),  and  (14)  constitute  a  set  of  four  source  functions 
in  a  counter- st reaming  beam-plasma  device.  These  four  sets  of  collisional 
interaction  mechanisms  are  best  illustrated  in  figure  6. 


Figure  6.  Four  sets  of  possible  enhance  Raman  scattering  interactions  in 
a  counter-streaming  beam-plasma  device. 


As  a  consequence,  there  will  be  four  sets  of  possible  enhanced  Raman 
scattering  interactions  in  a  counter-streaming  beam-plasma  device.  In  as  much 
as  these  four  source  functions  are  linearly  independent,  their  individual  con¬ 
tributions  to  the  Raman  scattered  transverse  electromagnetic  radiation  power 
output  from  the  device  are  also  independent,  and,  hence  they  may  be  calculated 
separately.  The  volume  emissivities  for  the  enhanced  coherent  Raman  power 
densities  Q(  C4jpe)  can  be  calculated  accordingly.  For  non- isothermal 
plasma,  Prasad  and  Lei  by  have  obtained  the  following  results: 


Q5  (~2  <V  Mo  3?  a  Fl(a>  “o  .  <'®> 


for  a  weakly  pumped  plasma,  and 

QT_(~2 _2JEe)  . 

Qr(~71v  "if" 


. ..  L.<  '111 ) 

35  a  Fj (a) 


H?  W2  +  W2W1 


(19) 


for  a  strongly  pumped  but  stable  plasma. 


is  the  power  density  of  a  single  beam-plasma  system. 


where 


is  the  total  power  density  of  a  counter-streaming  beam- 
plasma  device. 


H 


is  the  energy  density  of  the  pump-wave  in  the  single-beam 
system. 


W,  2  is  the  electrostatic  field  energy  density  for  the  two 
’  anti -parallel  pump  waves. 


F-j  (a) 


are  the  normalized  spatial  integrals,  with  ( a )/p 

between  2  and  3  for  physically  meaningful 
value  of  a. 


(a) 


a  =  Vc 


The  maximum  interaction  power  density  for  physically  realizable  counter¬ 


streaming  beam-plasma  device  is  given  by: 


QLx  (~2  ‘‘W  =  °*229(  3*  +  4a_1)  (mi  Te1}  ^  (L  Jb  Pb)3/2  . (20) 

for  the  case  of  Tg/T..  =  20,  V/<j=  0.1,  where  Tg  and  are  the  electron  and 
ion  temperatures  respectively.  Vis  the  electron  collision  frequency  and  nu  is 
the  mass  of  the  ions.  The  parameters,  oi  =  vph/vTe>  L,  Jb»  and  Pb  have  been 
defined  before. 
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III.  Objective  and  Design  Considerations 

Millimeter  wave  tube  employing  beam-plasma  interaction  mechanism  is  still 
in  its  infant  stage.  Little  work  has  been  done  both  theoretically  as  well  as 
experimentally.  There  are  several  problem  areas  that  must  be  tackled  and  solved 
before  a  feasible  engineering  model  can  be  developed.  Among  them  the  more 
immediate  concern  is  that  of  how  the  emitted  electromagnetic  radiation  can  be 
collected  and  re-direct  out  of  the  beam-plasma  interaction  volume  so  that 
frequency,  power,  or  other  measurements  can  be  made  external  to  the  beam- 
plasma  volume. 

Since  the  Raman  scattered  transverse  electromagnetic  waves  emitted  from 
the  beam-plasma  interaction  are  readily  radiated  out  of  the  plasma  volume, 
Prasad-Leibys '  devices  are  power  oscillators.  The  emitted  electromagnetic 
radiation  can  be  simply  collected  and  re-directed.  In  reference  (1),  a  coni¬ 
cal  metal  horn  has  been  proposed  as  a  simple  means  by  which  the  beam-plasma 
interaction  theory  for  the  generation  of  microwave  power  could  be  verified 
in  laboratory.  However,  in  many  engineering  applications,  this  may  not  be 
a  useful  proposition.  A  cavity,  or  waveguide  of  some  sort,  must  be  incorpor¬ 
ated  into  the  device  so  that  the  generated  microwave  power  can  be  controlled 
and  re-directed  to  meet  the  various  applications. 

It  is  the  objective  of  this  research  to  examine,  theoretically,  some 
othe  ways  and  means  that  may  be  more  viable  from  a  system  design  point  of  view. 
To  start  with,  let  us  consider  the  simpler  single  beam-plasma  model  first.  The 
emitted  electromagnetic  radiation  field  produced  by  this  model  is  that  of  the 
quadrupole  distribution  pattern  shown  in  fig.  3.  We  observe  that  the  forward 
radiation  pattern  is  very  sensitive  to  the  parameter  change  of  the  value  of  ct 
(=  v  h/c).  At  a  certain  value  of  o l  (=  0.6),  the  emitted  electromagnetic 
radiation  concentrates  in  a  very  narrow  conical  shape  that  is  almost  parallel 
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to  the  direction  of  the  electron  beam.  In  view  of  its  very  small  escape-angle, 
such  a  wave,  or  ray,  may  not  be  able  to  escape  from  the  plasma  column.  As  the 
value  of  oi  vthat  is,  the  beam  voltage)  gets  smaller  the  emitted  electromagnetic 
radiation  tilts  upward  and  its  escape-angle  becomes  larger.  This  will  enable  the 
ray  to  escape  from  the  plasma  column  depending  upon  the  dimensions  of  e  plasma 
column.  In  the  case  of  back  radiation,  the  situation  is  entirely  different. 

In  this  case,  in  view  of  their  relatively  large  escape-angle,  0^,  all  emitted 
electromagnetic  radiations  from  the  device  are  capable  of  escaping  from  the 
plasma  column  for  all  values  of  a  .  Furthermore,  these  escape-angles  show 
relatively  little  change  for  a  wide  variation  of  the  value  of  a  .  Thus,  if  « 
is  chosen  to  be  0.6,  corresponding  to  a  beam  voltage  of  92-kv.,  the  only 
emitted  electromagnetic  radiation  that  are  capable  of  escaping  from  the  plasma 
column  will  be  entirely  iR  the  backward  direction.  This  is  illustrated  in 
fig.  7. 


Figure  7.  Backward  radiation  pattern  of  a  single  beam-plasma  device 
having  oi  *  o.6. 
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With  reference  to  fig.  7,  from  eq.  (7),  the  electric  field  vector  Ey  is  in 
the  direction  of  the  vector,  n  x  k^,  where  fi  is  the  unit  vector  specifying  the 
direction  by  which  the  Raman  scattered  electromagnetic  radiation  is  to  be 
observed,  Jk^  is  the  direction  of  the  right-directed  electron  beam.  Thus,  the 
electric  field  lines  are  closed  circles  with  centers  located  along  the  axis  of 
the  cylindrical  electron  beam.  At  a  given  point  on  this  circle  the  electric 
field  vector,  which  is  tangent  to  the  circle,  is  pointed  in  the  anti-clockwise 
direction  when  viewed  along  the  direction  of  the  electron  beam.  Hence,  the 
electric  field  is  transverse  to  the  direction  of  the  electron  beam,  and  it  has 
the  axial  circularly-symmetrical  property.  Next,  by  eq.  (8),  the  magnetic  vector. 
By,  which  is  in  the  direction  of  the  vector  n  x  Ey,  is  perpendicular  to  both  the 
vectors,  n  and  Ey.  Thus,  at  points  P  and  Q  (fig.  7.)  the  vector  Ey  is  normal  to 
the  plane  of  the  paper,  whereas  the  vector  By,  which  is  normal  to  both  n  and  Ey, 
must  lie  also  in  the  plane  of  the  paper  at  these  two  points.  When  expressed  in 
cylindrical  polar  coordinates  the  magnetic  vector  By  may  be  decomposed  into 
parallel  and  radial  components.  The  parallel  components  Bz>  are  all  pointed  in 
the  negative  z-axis  direction,  whereas  the  radial  component,  B^,  are  all  directed 
radially  inward  toward  the  centers  of  the  electric  field  circles.  Therefore, 
the  emitted  electromagnetic  radiation  in  the  backward  direction  of  a  single 
beam-plasma  device  is  of  the  circularly-symmetrical  transverse  electric  type  - 
much  like  that  of  the  TEQ^  mode  in  a  circular  cylindrical  waveguide. 

For  o i  small  than  0.6,  the  escape-angle,  0^,  of  the  emitted  electromagnetic 
radiation  in  the  forward  direction  becomes  larger  as  the  escaped  rays  are  tilted 
upward.  As  <*  reduces  to  0.3,  corresponding  to  a  beam  voltage  of  23-kv.,  the 
emitted  electromagnetic  radiation  are  ready  to  escape  from  the  plasma  column 
if  the  length  and  the  size  of  the  plasma  column  were  suitably  chosen.  By  a 
similar  field  analysis,  it  can  be  shown  that  the  forward  radiation  field  pattern 


also  possesses  axial  circularly-symmetrical  transverse  electric  property.  The 
polarity  of  the  circular  electric  field  lines  is  the  same  as  that  of  the  back 
radiation  field.  Although  the  parallel  components  of  the  magnetic  field,  Bz, 
has  the  same  polarity  as  that  which  produced  by  the  back  radiation  magnetic 
field,  but  its  radial  component  field,  Br,  has  reversed  in  phase  by  180  degrees. 
This  is  illustrated  in  fig.  8. 


Figure  8.  Forward  and  backward  radiation  patterns  of  a  single  beam-plasma 
device  having  &  =  0.3. 

No  precise  field  data  have  been  given  by  Prasad  and  Leiby  concerning  their 
counter-streaming  device,  but  in  light  of  what  has  been  discussed  in  section  II, 
the  Interaction  mechanism  of  this  device  may  produce  a  much  more  complicated 
field  pattern  than  that  shown  in  fig.  3.  However,  it  is  conceivable  that  the 
field  distribution  of  a  counter-streaming  device  may  still  be  of  the  circularly- 


symmetrical  transverse  electric  type  capable  of  emitted  electromagnetic  radiation 
in  both  the  forward  and  backward  directions  -  much  like  that  of  fig.  8. 


This  field  characterization  provides  us  much  insight  information  as  t.o 
how  the  microwave  collector  system  should  be  designed.  There  are  other  problem 
areas  to  be  considered  too.  In  a  beam-plasma  device,  large  interaction  volume 
is  required  for  the  realization  of  the  very  large  power  generating  capability 
of  the  device.  One  is  led  to  consider  a  cylindrical  cavity  or  waveguide  that 
is  large  compared  to  the  wavelength  contemplated.  The  oversize  cavity  or  wave¬ 
guide  may  support  a  large  number  of  unwanted  modes  within  the  frequency  range 
of  interest.  Beam-plasm a  device  using  multi-mode  cavity  or  waveguide,  may 
encounter  considerable  difficulty.  For  instance,  oscillations  may  be  set  up 
in  one  mode,  then  the  others,  or  a  continual  change  of  modes,  resulting  in  a 
sudden  jump  in  frequency.  Furthermore,  since  the  density  distribution  of  the 
plasma  can  not  remain  undisturbed  in  the  presence  of  a  strong  cavity  or  wave¬ 
guide  fields  within  the  interaction  plasma  volume,  a  very  small  change  in 
density  distribution  may  produce  a  shift  in  oscillations  from  one  mode  to 
another  resulting  in  a  considerable  variation  in  frequency.  Therefore,  the 
plasma  tube  should  not  be  placed  in  the  high  field  region  of  either  the  cavity 
or  waveguide,  whichever  microwave  collector  system  is  used. 

Another  problem  associated  with  a  beam-plasma  device  using  Raman  scattering 
interaction  mechanism  is  the  possibility  of  the  emitted  electromagnetic  radiation 
re-interacting  with  the  plasma  upon  reflection  from  the  cylindrical  metal,  or 
dielectric,  wall  of  the  cavity  or  waveguide.  Since  the  emitted  electromagnetic 
radiation  frequency  2  &Jpe)  is  much  greater  than  the  plasma  frequency,  <Ope, 

the  plasma  is  transparent  to  the  return  ray,  and  the  re- interaction  of  the  return 

ray  with  the  plasma  may  produce  adverse  effect  which  may  impair  the  proper 

operation  of  the  device.  Hence  it  should  be  avoided  by  all  means. 


_ 


_ . 


-  25  - 


Based  on  the  afore  mentioned  guide-lines,  for  a  proper  operation  of  a  beam- 
plasma  device,  three  microwave  collector  structures  have  been  proposed  and  will 
be  investigated  in  detail  in  the  subsequent  sections.  They  are: 

Type  I:  The  hollow  dielectric  tube  waveguide. 

Type  II:  The  cylindrical  metal  waveguide  partially  filled  with  a  coaxial 
dielectric  tube. 

Type  III:  The  hollow  dielectric  tube  waveguide  supported  by  a  conducting 
cylinder. 

All  three  structures  can  be  considered  as  the  variant  of  the  basic  dielectric 
tube  transmission  line.  In  its  most  elementary  form,  it  consists  of  a  dielectric 
tube  immersed  in  air  as  shown  schematically  in  fig.  9.  The  dielectric  tube  (of 
inner  and  outer  radii  "a"  and  "b"  respectively)  separates  the  surrounding  space 
into  three  regions,  denoted  by  I,  II,  and  III,  respectively.  We  denote  the 
dielectric  permitivities  of  the  three  regions  by  6^,  62>  and  respectively, 
and  assume  also  that 

Region  III 


Pig. 9.  Schematic  diagram  of  a  dielectric  tube 
waveguide  structure. 


1 )  Region  I :  0  <  r  <  a 

This  is  the  interior  of  the  dielectric  tube  where  the  cylindrical  beam- 
plasma  tube  proper  is  to  be  housed.  The  plasma  tube  may  be  tight-fitted 


or  loosely-fitted  into  the  hollow  empty  space.  For  the  loosely-fitted 
case,  a  separate  dielectric-tube  sleeving  of  small  dielectric  constant 
may  be  used  to  hold  the  plasma  tube  in  place  so  as  to  attain  the  desired 
axial  symmetrical  property  of  the  complete  unit. 

2)  Region  II :  a  c  r  <  b 

This  is  the  region  of  space  whereby  the  generated  microwave  energy  of 
the  device  is  to  be  contained,  collected,  and  re-directed.  This  would 
require  the  dielectric  material  used  to  be  of  the  very  lossless  type 
with  a  dielectric  constant  larger  than  that  used  in  both  the  regions 
I  and  III. 

3)  Region  III:  r  >  b 

This  may  be  simply  the  free  space  as  in  the  case  of  the  type  I  structure, 
or  it  may  be  a  grounded  conducting  cylindrical  tube  as  in  the  case  of 
type  II  structure. 

Type  II  structure  may,  therefore,  be  considered  as  a  modified  form  of  the 
type  I  structure  in  which  region  III  is  replaced  by  a  grounded  conducting  cylin¬ 
drical  tube  having  an  inner  radius  "b". 

Type  III  structure  may  also  be  considered  as  a  special  case  of  the  type  I 
structure  in  which  region  I  is  being  replaced  by  a  grounded  conducting  cylinder 
having  an  outer  radius  "a".  All  three  microwave  collector  structures  are  to  be 
operated  in  the  circularly-symmetrical  transverse  electric  modes.  Field  analysis 
indicates  that  the  first  two  structures  possess  several  very  desirable  properties 
in  millimeter  wave  applications.  First,  these  two  structures  can  be  readily 
integrated  into  the  beam-plasma  tube  proper  to  form  a  complete  unit  like  most  of 
the  conventional  microwave  tubes,  such  as  the  traveling-wave  tubes.  Secondly, 
these  two  structures  are  quite  capable  of  confining  the  emitted  electromagnetic 
radiation  fields  within  Its  high  dielectric  constant  medium.  Thirdly,  both  these 
structures  are  easily  adopted  to  function  as  microwave  collector  as  well  as  to 
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re-direct  the  microwave  energy  out  of  the  beam-plasma  unit  when  one  end  of  the 
structure  is  closed  by  a  grounded  metal  plate  and  the  other  end,  which  may  be 
appropriately  tapered  to  a  suitable  size,  is  left  open  to  let  the  microwave 
energy  leaving  the  device.  Obviously,  by  the  lacking  of  these  confining, 
collecting,  and  re-directing,  capabilities,  type  III  structure  is  not  as 
versatile  as  the  other  two  structures.  Nevertheless,  type  III  structure  does 
have  limited  applications  as  a  millimeter  wave  transmission  line. 

The  theoretical  foundation  and  design  parameters  of  these  microwave 
structures  will  be  presented  in  sections  IV,  V,  and  VI.  In  section  IV,  the 
general  formulation  of  dielectric-tube  transmission  line  will  be  given  first. 
Solution  of  type  I  transmission  line  will  be  followed  with  results  analyzed 
and  discussed.  Section  V  deals  entirely  with  the  solution  of  the  type  II 
transmission  line  to  be  followed  by  a  detailed  discussion.  Since  type  III 
transmission  line  has  only  a  limit  application,  its  solution  is  given  briefly 
in  section  VI. 

IV.  The  Hollow  Dielectric  Tube  Waveguide. 

4a.  Introduction 

Dielectric  waveguides  are  here  defined  as  primary  transmission  lines 
which  employ  a  system  of  dielectric  elements,  as  distinct  from  a  system  of 
conductors,  such  as  coaxial  transmission  lines  or  metal  waveguides,  to  transmit 
electromagnetic  energy  from  one  or  more  directions.  In  such  waveguides 
electromagnetic  waves  are  guided  along  the  axis  of  a  cylinder  at  the  surface 
of  which  exists  a  discontinuity  in  the  dielectric  properties  of  the  mediums 
inside  and  outside  the  cylinder.  A  rod  of  dielectric  material  is  an  example 
of  such  a  waveguide,  as  is  also  a  cylindrical  hole  in  a  block  of  dielectric 
material. 
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The  propagation  of  electromagnetic  waves  along  a  solid  dielectric  cylinder 
was  investigated  theoretically  by  Hondros  and  Debye^®^  early  in  1910.  Some 
experimental  work  has  been  done  by  Zahn^19^  and  Schriever^20^  during  the  period 
from  1916-1920.  Later  in  1936,  the  subject  was  further  advanced  by  Barrow^21  \ 
Southwork^22^ ,  Carson,  Mead,  and  Schelkunoff ^23^ .  During  the  World  War  II, 
the  development  of  the  metal  tube  waveguides  made  rapid  progress  and  as  trans¬ 
mission  systems  completely  overshadowed  the  dielectric  waveguides,  largely  due 
to  the  fact  that  the  electromagnetic  field  is  contained  entirely  within  the 
matal  tube  guide  while  it  also  exists  outside  the  dielectric  guide  and  leads 
to  greater  transmission  loss  due  to  radiation  when  bends  and  discontinuities 
were  present  in  the  transmission  line.  During  the  period  from  1940-1950, 
most  of  the  work  done  was  in  the  field  of  dielectric  rod  antennas.  The  works 
of  Southwork^24a^ ,  Mueller  and  Tyrell^24^,  Halliday  and  Kiely^24c\  Watson 
and  Horton (24dM24e),  simon^24^,  Wilkes^24^,  Horton,  Karal ,  and  McKinney^24^ 
and  others,  have  brought  the  theory  of  the  subject  of  radiation  from  dielectric 
rods  to  its  most  advanced  state.  Published  work  on  dielectric  tube  antenna  is 
much  less  extensive  than  that  on  dielectric  rod. 

The  problem  of  guided  wave  propagation  along  a  dielectric  tube  has  been 

considered  theoretically  by  Zachoval^2^,  in  1932.  In  1934  Liska^33^  published 

the  results  of  his  experiments  made  to  confirm  the  theory  of  Zachoval .  He 

considered  on  the  TM~  modes  but  used  several  values  of  wall  thickness.  His 

Op 

results  were  in  good  agreement  with  Zachoval's  theory. 

The  most  complete  treatment  of  the  theory  of  the  dielectric  tube  wave¬ 
guide  was  carried  out  in  1949  by  Astrahan^34^ ,  who  solved  the  characteristic 
equations  for  TMq^ ,  TEq^ ,  and  the  hybrid  HE^  modes.  He  obtained  curves  of 
the  variation  of  wavelength  of  the  guide  wave  with  tube  dimensions  and  frequency 
of  the  exciting  radiation.  His  theoretical  results  were  substantiated  by  experi 
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mental  results.  However,  Astrahan's  systematic  trial  and  error  method  of 
solution  of  the  characteristic  equation  is  rather  unwieldy  and  tedious.  We 
shall  seek  the  modern  computer's  help  to  obtain  from  the  solution  of  the  char¬ 
acteristic  equation  some  relevant  informations  and  data  which  have  direct 
implication  in  our  millimeter  wave  applications. 

In  this  section  we  will  pay  our  main  attention  to  the  principle  character¬ 
istics  of  the  guided  wave  propagation  in  a  dielectric  tube.  The  dependence  of 
the  waves  on  dielectric  tube  dimensions  and  dielectric  constant  for  the  circularly- 
symmetrical  transverse  electric  modes  will  be  derived  theoretically  and  those 
features  of  wave  propagation  along  the  dielectric  tube  which  have  a  bearing  on 
the  design  of  dielectric  tube  microwave  collector  will  be  pointed  out.  The 
circularly-symmetrical  transverse  magnetic  and  the  asymmetric  modes  will  not 
be  dealt  with  by  reasons  discussed  in  section  III. 

4b.  Wave  Propagation  Along  a  Dielectric  Tube. 

Consider  a  cylindrical  dielectric  tube  having  inner  and  outer  radii 
"a"  and  "b"  respectively.  The  dielectric  tube  separates  the  space  surrounding 
it  into  three  regions  having  different  dielectric  permitivities,  £'s,  and 
magnetic  permeabilities,  jj'  s.  It  is  also  assumed  that  these  dielectric 
mediums  are  isotropic  such  that: 

(a)  -  -^3  *  where  jJQ  is  the  absolute  magnetic  permeability 
of  the  vacuum, 

(b)  6Z  >  63  =  €-| ,  and 

(c)  real  conductivity  in  all  three  mediums  is  zero. 

The  cylindrical  co-ordinate  system  used  in  the  subsequent  analysis  is  shown 
In  fig.  10,  where  the  co-ordinates  of  a  point  P  in  the  cylindrical  space  are 
(r,  0,  z). 


Fig. 10.  Cylindrical 

co-ordinate  system 


It  is  further  assumed  that: 

1)  Rationalized  M.K.S.  system  of  units  is  used. 

2)  Waves  are  guided  along  the  positive  z-axis  direction. 

3)  The  z  and  t  dependence  of  all  field  quantities  is  written  as  e( jco t  -flz), 

where  u  is  the  radian  frequency  of  the  wave  and  is  the  propagation 
constant  which  may  be  complex  or  pure  imaginary.  With  this  dependence, 

we  can  replace  d/  qz>  and  3/0t  in  Maxwell's  equations  by  -  $  and  jc») 
respectively. 

For  charge-free  mediums  with  zero  conductivity.  Maxwell's  equations  can  be 
written  as 

V  x  E  =  -  jo)  H  (4-la) 

y  xH  =  j  a)  €  E  (4-lb) 

*7  *  E  =  0  (4-lc) 

V  •  H  =  0  (4-1 d ) 

A  time  dependence  of  j  a)  t  is  assumed. 

Defining  the  unit  vector  az  as  the  direction  of  wave  propagation,  then,  we 
can  express  all  field  components  transverse  to  uz  in  terms  of  the  longitudinal 
field  components  parallel  to  Thus,  we  let 
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I  =SzE2+It 

H  =  a2  Hz  +  Ht 


(4-2) 


where  Ey  and  Hy  are  the  transverse  components  of  the  electric  and  magnetic  field 


vectors  respectively.  We  obtain  from  the  first  curl  equation,  (4-la) 
(  VT  -  az)  x  (  Et  +  Ez  az)  =  -  j*>  PQ  (  Hy  +  Hz  ^) 

where  the  "del"  vector  operator  V  has  been  written  as 

V=7r*  Aa  -i- 
d  z 

and  the  transverse  "del"  operator  Vy  is  given  by 

VT  =  a  ^  JL  A_ 

T  r  3  r  *  r  d* 


(4-3) 


(4-4a) 


(4-4b) 


in  cylindrical  co-ordinates. 


Eq.  (4-3)  can  be  separated  into  transverse  and  longitudinal  component 
vector  equations  as  follows: 

a)  The  transverse  component  vector  equation: 

Vy  x  Ez  az  -  3  az  x  Ey  =  -  jcd  Hy  (4-5a) 

b)  The  longitudinal  component  vector  equation: 

VT  x  Ey  =  -  j  co  Hz  az  (4-5b) 

In  like  manner,  in  the  absence  of  real  conductivity,  the  second  curl 
equation  (4-lb)  can  also  be  separated  into  transverse  and  longitudinal  compo¬ 
nent  vector  equations  as  follows: 

^7y  x  Hz  az  -  az  x  Hy  =  j6)£  Ey  (4-6a) 

and 

Vy  x  Hy  =  j  c»>  €  Ez  az  (4-6b) 

By  an  appropriate  manipulation  of  vector  algebra  we  can  express  the  two 
transverse  components,  Hy  and  Ey,  in  terms  of  the  two  longitudinal  components, 

Hz  and  Hz»  This  is  done  by,  first,  substituting  Ey  from  eq.  (4-6a)  into  eq.  (4-5a) 


' 


to  yield 


Ht  =  jo)  6 
T  h2 


( VT  x  E2  a2)  -  <  VT  H. 

n 


(4-7) 


o  p  p  o  p 

where  h  =  k  +  "3  and  k  =  6) JUQ  6  ,  is  the  dispersion  relation  for 
isotropic  mediums. 

Next,  substituting  eq.  (4-7)  into  eq.  (4-6a)  we  obtain 


%  ■  <  V  Hz  4,)  -  i  VT  Ez  (4-8) 

where  h4  is  the  same  as  that  which  is  given  in  eq.  (4-7). 

Thus,  knowing  the  two  longitudinal  component  fields,  Hz  and  Ez,  the  two 

transverse  component  vectors,  and  Ey,  can  be  readily  obtained  from  eqs.  (4-7) 
and(4-8)  respectively.  Therefore,  the  whole  field  problem  is  now  resolved  into 
the  question  how  the  two  longitudinal  component  fields,  Hz  and  Ez>  can  be  deter¬ 
mined?  This  can  be  done  by,  first,  substituting  eq.  (4-7)  into  eq.  (4-6b),  and 
after  appropriate  vector  manipulation  we  obtain 

(  Vy  +  h2  )  Hz  =  0  (4-9) 

2 

where  \7y  is  the  transverse  Laplace  operator.  And  in  like  manner,  substituting 
eq.  (4-8)  into  eq.  (4-5b)  we  obtain 


(  Vj  +  h2  )  Ez  =  0 


(4-10) 


Eqs.  (4-9)  and  (4-10)  are  homogeneous  Helmholtz  equations  for  Hz  and  Ez  respectively. 

p 

By  expressing  the  transverse  Laplacian  operator  VT  in  cylindrical  form,  eqs.  (4-9) 
and  (4-10)  become 


( 1 

\  r  dr 


(  r  d_  )  +  1 

dr 


-  *■  r 

3-  +  h‘ 

a  u2 


(4-11) 


2  2  2  2  2 

with  h  =  k  +  3  »  where  k  -  OJ  /JQG  •  When  the  longitudinal  component  fields, 
and  Ez»  are  solved  from  eqs.  (4-11)  the  transverse  component  fields,  and  Ej, 
can  be  determined  from  eqs.  (4-7)  and  (4-8)  respectively. 

Eqs.  (4-11)  may  be  solved  by  the  standard  method  of  separation  of  variables. 


For  example,  to  solve  eq.  (4-11)  for  Hz  we  let 


H.(r,#)  *  R(r)  #  (0) 


(4-12) 
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where  R(r)  is  a  function  of  only  r,  and  <$(0)  is  a  function  only  of  0.  Then, 
by  performing  the  needed  partial  differentiation  we  obtain  the  following 
equations: 


d2$  =  -/\2f  (4-13a) 

d7 


and 


+  ldRf(h2-XL)  R  =  0 
r  dr  1 


(4-1 3b ) 


2 

where  A  is  a  separation  constant  which  may  be  any  real,  imaginary,  or  complex 
number. 

The  solution  of  eq.  (4-13a)  is  sinusoidals,  which  can  be  the  combination 
of  sin7\0,  cosa0,  or  e±J  however,  in  view  of  the  fact  that  fields  inside 
a  cylindrical  dielectric  waveguide  having  a  circular  cross-section  must  be  a 
periodic  function  of  0  with  period  in  2 7T ,  therefore,  the  constant  A  must  be 
an  integer.  We  shall  choose  *  =  n.  Replacing  *  by  n  in  eq.  (4-13b),  and  by 
a  transformation  of  variables  eq.  ( 4-1 3b)  becomes 


A 

-r 


+  1  dR  +  (1  -  jl  )  R  =  0 

1  ^  r* 


(4-l3c) 


where  'f  =  hr,  and  h  is  defined  as  before.  Eq.  ( 4- 1 3c )  is  the  Bessel  equation 
of  order  n.  The  most  general  solutions  of  eq.  (4-1 3c)  are  linear  combinations 
of<25> 

R(r)  *  F  Zn(  f  ) 

where  F  is  an  arbitrary  constant,  is  the  cylindrical  functions  of  order  n, 

which  may  be  Bessel  function  Jn('f),  Neumann  function  Nn(^),  Hankel  function  of 
the  first  kind  (f)»  or  Hankel  function  of  the  second  kind  Choice 

of  a  particular  function  will  depend  upon  the  region  of  interest  considered.  We 
will  discuss  it  in  detail  in  the  next  section. 

With  both  $  and  R  solutions  obtained,  the  product  of  R  and  $  gives  us  the 


_ -  , 
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desired  solution  for  the  longitudinal  magnetic  field 
H2(r.0)  •  F  Zn(hr)  <“*  *) 

In  like  manner,  we  can  solve  eq.  (4-1 0)  co  yield 


(4-14) 


(4-15) 


Ez(r,0)  =  G  Zn  (hr)  *> 

where  G  is  another  constant.  The  propagation  factor  e  z  is  not  included  in 

both  eqs.  (4-14)  and  (4-15). 

Finally,  with  both  H and  Ez  determined,  the  transverse  electric  and  magnetic 
fields  may  be  obtained  from  eqs.  (4-8)  and  (4-7).  The  result  is  given  in  the 
following: 

E.  (r,t)  -  -  I'M  (  l  4^  )  - 


JL  ^ 
h2  Tr' 


E0  (r ,0)  =  jft)^ 

Jl  d  r 


1 

r 


)_[z 

h 


(4-16a) 
( 4-1 6b ) 


and 


Hr  (r,0)  =  j  6)  €  (  l  dL)  -  Jf_  AHz 

h2  r  h2  <)r 


H0  (r,0)  =  -  jtJ  € 

h2 


JL  (1  AH1) 

h2  r  dt 


(4-1 7a ) 
(4-1 7B) 


This  completes  the  general  solution  of  the  wave  propagation  along  a  cylin¬ 
drical  dielectric  tube  waveguide  problem.  However,  based  on  the  guided  wave 
theory,  several  guide-wave  modes  are  possible  for  a  given  guiding  structure 
depending  upon  the  boundary  conditions  imposed.  For  instance: 

1)  The  TE  or  the  H-modes,  in  which  all  electric  fields  are  transverse  to 
the  direction  of  the  wave  propagation.  In  the  absence  of  the  longi¬ 
tudinal  electric  field,  Ez  *  0,  field  solutions  that  satisfy  all 
boundary  conditions  Imposed  by  the  guiding  structure  are  derived  from 
the  longitudinal  magnetic  field,  Hz,  alone.  This  is  readily  obtained 


- I — I — - I - u - m  I _  -  . I . I  -  -  ■  _ — _ . . _ I 
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from  the  solution  of  eq.  (4-9)  to  yield  Hz  as  given  in  eq.  (4-14). 

2)  The  TM  or  the  E-modes,  in  which  all  magnetic  fields  are  transverse  to 
the  direction  of  the  wave  propagation.  With  H z  =  0,  field  solutions 
that  satisfy  all  bcndary  conditions  imposed  by  the  guiging  structure 

are  derived  from  Ez  alone.  This  is  readily  obtained  by  solving  eq.  (4-10) 
to  yield  Ez  as  given  by  eq.  (4-15). 

3)  The  HE  or  the  hybrid  modes,  in  which  both  Hz  and  E?  are  required  to 
satisfy  Maxwell's  equations  and  the  boundary  conditions. 

4)  The  TEM  modes,  in  which  all  electric  and  magnetic  fields  are  transverse 
to  the  direction  of  propagation  with  both  Ez  and  Hz  are  zero.  In  these 
modes  both  the  transverse  fields,  E-j-  and  can  be  derived  from  the 
potential  functions  that  satisfy  the  Laplace  equation  instead  of  the 
Helmholtz  equation.  For,  if  ^z  =  Hz  =  °»  from  ec)s-  ( 4-5b)  and  (4-6b), 
we  obtain 

V  -j-  *  =  0  and  Vy  x  My  =  0 

Therefore,  both  Ey  and  are  irrotational ,  and  by  Helmholz  theorem, 

ET  =  -VT$e  and  Ht.-Vt#„ 

where  and  are  the  potential  function  for  the  electric  and  magnetic  fields 
respectively.  Furthermore,  in  a  source-free  region,  Vy  ’  iy  =  0  and 
Vr  ‘  H,  =  0.  This  gives  us  the  required  Laplace  equations: 

Vt  #3  and  Vt  #„  ‘°- 

As  discussed  in  section  III,  the  Prasad-Leibys 1  devices  generate  a  circularly- 
symmetrical  transverse  electric  waves.  The  boundary  conditions  of  all  three  micro- 
wave  structures  proposed  in  section  III  can  be  met  by  the  TE  or  H-modes  solution 
described  above.  In  what  follows  and  in  subsequent  sections  V  and  VI,  we  shall 
concern  ourselves  only  with  this  type  of  wave  solution  also. 
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4c.  Circularly-Symmetrical  Transverse  Electric  Field  Distribution. 

For  circularly-symmetncal  transverse  electric  modes  which  are  independent 
of  the  angle  variable  0,  the  field  components  of  this  type  of  waves  travelling  in 
the  positive  z-direction  can  be  derived  from  the  longitudinal  magnetic  field  Hz, 
as  discussed  in  section  4b.  The  solution  of  the  Helmholz  equation,  (4-9),  for  Hz 
yields  the  desired  field  equation,  (4-14).  Let  us  consider  the  most  general  form 
of  eq.  (4-14)  for  each  of  the  three  regions. 

In  regions  I  and  II,  the  cylindrical  functions  Zn  (hr)  may  be  written  as 
the  combinations  of  the  Bessel  functions,  Jn  (hr)  and  Neumann  functions,  Np  (hr). 
However,  in  region  I  (0<r  <a)  which  includes  the  origin  of  the  cylindrical  co¬ 
ordinates,  the  Neumann  function  can  be  expressed  in  the  form 


Nn  (hr) 


r  (n)  (  2  ) 


rr 


h  r 


en 


0 


In  order  to  ensure  a  finite  field  at  r  =  0,  Nn  (hr) 


which  is  infinite  as  r — *0 
functions  must  be  excluded  from  the  solution.  Then,  for  region  I,  0  r<  a, 
the  following  is  a  suitable  form  for  eq.  (4-14). 


Hzl(r,0,z)  -F  J,(h,r)  (i“J5)  e-  * z  (4-18) 

io)  t 

The  time  factor  of  e°  has  been  assumed. 

In  region  II,  a  ^  r  <  b,  which  does  not  include  the  origin,  both  the  Bessel 
and  Neumann  functions  must  be  included,  the  general  solution,  eq.  (4-14),  takes 
on  the  form: 


Hz2(r,tf,z)  =  fA  Tn 


(h2r) 


Nn(h2r) 


(cos 
sin 


n0 

n0 


Yz 


(4-19) 


In  region  III,  r  >  b.  This  is  the  region  outside  the  dielectric  tube,  the 
appropriate  cylindrical  function  to  use  is  the  combinations  of  the  Hankel  functions 


(h,r)  and  (h.r),  where 
n  j  n  j 

hJ”  <z)  =  Jn  (z)  +  i  Nn(z) 
»n2>  <2>  ‘  Jn  (4)  *  J  N„(z) 
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are  the  Bessel  function  of  the  third  kind,  or  the  Hankel  function.  The  choice 
is  based  on  the  following  reasoning.  In  region  III  the  dispersion  equation  is 


written  as 


=  k*  ♦  f 
=  (  G>  )2 


(  JuL  ) 
ug 


(4-20) 


For  propagation  mode,  f  -  jp  .  p  =  &tyug,  where  u^  is  the  velocity  of 
propagation  of  the  guided  wave  along  the  dielectric  tube,  and  k3  =  cj  (  jjq  £3)* 
and  u  ,  =  (  p  €,)"*,  is  the  velocity  of  uniform  plane  wave  in  medium  III.  Since 

we  have  assumed  €2  •-*  €3,  up2  (  =  (  juq  £2)~*  )  is  1ess  than  up3;  but  ug»  the 

velocity  of  the  guided  wave,  or  the  group  velocity  of  the  guided  wave,  can  not  be 

9 

greater  than  its  phase  velocity.  Eq.  (4-20)  implies  that  h3  must  be  less  than  or 
equal  to  zero,  and  so  h3  must  be  either  zero  or  pure  imaginary. 

Let  us  examine  further  the  asymptotic  behavior  of  these  Hankel  functions  for 


large  r.  For  large  r _ 

H  0)  fhrl /3Z  ej(hr  “  J2JL  -  _IL  ) 

Hn  r  e  2  ~ 


Hn(2)  (hr)- 


- j  (hr  -  n  T7  -  JT  ) 


(4-21) 


We  see  that  as  r->00,  behaves  as  e^r  and  H  as  e'^hr.  In  order  for 

waves  to  be  guided  inside  the  dielectric  tube  (a  <1  r  <  b)  the  fields,  which  must 
approach  zero  at  large  distance  from  the  axis,  must  be  envanescent  in  the  radial 
direction.  The  proper  choice  of  the  solution  for  region  III  will  be  Hankel 
function  with  an  imaginary  argument  known  as  the  modified  Hankel  function.  Such 
a  choice  will  ensure  the  proper  decrease  of  the  field  at  infinite  and  will  also 
give  the  proper  direction  of  power  flow  if  the  complex  parameter  h3  is  in  the 
first  quadrant  of  the  complex  plane.  The  most  general  combinations  of  Bessel 
functions  which  approach  zero  for  large  pure  imaginary  positive  argument  is 
Hn^  (hjr).  Hn^  (b3r)  must  be  excluded  from  the  solution  for  this  region. 


»  • 


Hence,  the  proper  formof  eq.  14-14  for  region  III  will  be 


Hz3(r,^,z)  =  C  Hn' 


<h3r) 


( cos  n0  ^ 
^sin  n0 ) 


(4-22) 


if  we  restrict  to  be  a  positive  pure  imaginary  quantity,  where  C  is  a  constant, 
and  h3  =  j  kc3,  with  kc3  a  positive  real  number. 

As  mentioned  earlier  we  are  seeking  the  ci rcularly-symmetrical  transverse 
electric  modes  solution  which  are  independent  of  the  angle  variable  0.  This  can 
be  satisfied  by  choosing  n  =  0.  Thus,  by  letting  n  =  0,  we  obtain  from  eqs.  (4-18), 
(4-19),  and  (4-22),  the  longitudinal  magnetic  field  in  the  three  regions.  They 
are  summed  up  in  the  following: 

Hzl  -  F  J0(h,r)  e-  *z 


Hz2  »  (A  J0(h2r)  +  B  N0(h2r)} 

H23  =  CHo(”  <h3r>  e'tZ 


1  CUT 

where  h3  =  j  kc3>  with  kc3  a  real  positive  number.  The  time  factor  eJ  has 


been  assumed. 


With  the  longitudinal  magnetic  field  of  the  three  regions  determined,  the 
transverse  electric  and  magnetic  fields  in  these  three  regions  can  now  be 
obtained  from  eqs.  (4-16)  and  (4-17).  The  result  is  listed  in  the  following: 


Region  I:  0<£  r  ^  a 
Hzi  =  F  J0(hlr)  e- 


H*1  =0 


F  J  •  (hlr) 


(24a) 


(24b) 


Erl  =0 


E01  =  F  Jq  (h-j  r)  e' 

”l 

r  _  n  1 


(4-24c) 


*f  » 
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Region  II:  a  <  r  <  b 

Hz2  =  (A  Jo(h2r)  *  8  Vh2r)}  e"^Z  (25-a) 

H*2  -  0 

Hr2  =  "  [  A  Jo(h2r)  +  8  No(h2r)  }  e'  ^  (25-b) 

Er2  =  0 

E*2  =  “R-”?  j  A  J0(h2r)  +  8  No(h2r)|  e'^Z  (25-c) 


where 


Region  III:  r  > 

b 

Hz3  =  C  Ho‘ 

(h3r)  e"^z 

(4-26a) 

v>  =  ° 

Hr3  ’  -  c  * 
h3 

.  H0('>  (h3r)  e  -®z 

(4-26b) 

Er3  =  0 

E,,  =  C  JCJ^o 

*3  -17- 

E23  =  ° 

H0(,)  (h3r)  e-*z 

(4-26c) 

h^  =  k^  +  3  2 

-  fr,  4  -  P2 

(4-27a) 

h^  »  kp  +  U2 

h3  =  kc3^ 

■  er2ko  -P2 

(4-27b) 

"  k3  +  ^ 

■  *r3ko  -P2 

(4-27c) 

where 

a)  %  is  identically  the  same  in  h^ ,  and  h^»  and  for  propagation  mode, 
^  *  !  P  »  where  /5  is  the  phase  shift  constant; 

b)  k  =  cJ  (  ^  6 o  £r)*  =  kQ  (  Gr)*t  and  kQ  =  cJ/c,  c  is  the  velocity  of 
light  and  £  p  Is  the  dielectric  constant; 

c)  the  prime  denotes  differentiation  of  the  function  with  respect  to  its 
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(25) 

argument,  and  the  recurrency  formula'  ’  for  the  cylindrical  functions, 

JQ(hr),  NQ(hr),  or  HQ^  (hr)  is 

z!  (hr)  =  d  Zo(hr)  =  -  Z, (hr)  (4-28) 

0  d(hr)  1 

and 

d)  F,  A,  B,  and  C,  are  constants  which  must  be  determined  by  the  boundary 
conditions  at  the  surfaces  r  =  a  and  r  =  b  such  that  the  tangential 
components  of  the  electric  and  magnetic  fields  must  be  continuous  at 
these  surfaces.  Thus 

(a)  Continuity  of  the  z-components  of  the  magnetic  fields  at  the  surfaces 
r  =  a  and  r  =  b  requires 

(i)  Hz1  =  Hz2  at  r  =  a,  and 

(ii)  Hz2  =  Hz3  at  r  =  b. 

This  will  give  us  from  eqs.  (24a),  (25a),  and  (26a) 

F  J0(h-| a )  =  A  J0(h2a)  +  B  NQ(h2a)  (4-29) 

C  H(1}(h3b)=A  Jo(h2b)  +  B  NQ(h2b)  (4-29') 

(b)  Continuity  of  the  tangential  components  of  the  electric  fields  at  the 
surfaces  r  =  a  and  r  =  b  requires 

(i)  E^-j  =  E^2  at  r  =  a,  and 

(ii)  E*3  =  E*2  at  r  =  b* 

These  will  give  us  from  eqs.  (4-24c),  (4-25c),  and  (4-26c) 

F  h'1  j’t^a)  =  h'1  [A  j‘(h2a)  +  B  (iy)]  (4.30) 

I 

C  h’1  H<])  (h3b)  =  h2][A  J^(h2b)  +  B  N;  (h2b)]  (4-30’) 

We  have,  thus,  obtained  an  array  of  four  equations  in  four  known.  In 
order  to  solve  the  problem  completely,  the  values  of  the  four  arbitrary  con¬ 
stants,  A,  B,  C,  and  F  must  be  determined.  The  expressions,  eqs.  (4-27)  to 
(4-30),  form  a  set  of  four  homogeneous  equations  in  these  unknown  coefficients. 
It  is  an  elementary  result  of  the  theory  of  linear  algebraic  equations  that 
such  a  set  of  equations  can  hay  no  solution  other  than  the  trivial  one 
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A=B=C=F=0  unless  the  determinant  of  the  coefficients  of  the  system 
vanishes.  This  implies  that  there  are  no  solutions  of  physical  significance 
unless  the  determinant  formed  by  the  coefficients  of  these  four  unknown  vanishes. 
Thus,  by  rearranging  eqs.  (4-27),  (4-28),  (4-29)  and  (4-30)  in  an  appropriate 
order,  we  have 


0o(h2a) 

0„(b2b) 

Jo(h2a> 

N0Ch2a) 

No(h2b> 

N0(h2a) 

0 

-Ho0)<h3b> 

0 

-  J0(V> 

0 

■(t)  J°(hia) 

=  0 

Jo<h2b> 

NJh2b) 

-(  ^)H‘"(h3b) 

\  hi  / 

0 

Expansion  of  this  determinant  leads,  if  we  discard  trivial  solution,  to  the 
following  transcendental  equation: 
uO) 


(jy)  0 


(jy) 


[ 


h0^7 


jn(x)  Nn(mx)  -  Jn(mx)Nri(x)  +  ($) 


J0(jmyh1h3"1 ) 


Jn(jmyh1h  ) 


|"jn(x)N'(mx)  -  J  (mx)N  (x'  1  +  (£) 

L°  0  0  0  i  y  J0{Jmyh1h3-« ) 

where  for  conciseness,  we  have  adopted  the  notations 
m  -  a/b,  x  •  h2b,  y  •  kc3b, 
with  h3  *  jkc3»  and 


kb 


*r2  ko 

(jk„J2 


5  f 


c3J 
k2 
r3  Ko 


-  ft* 


^j0(x)No(mx)  -  J0(mx)N0(xjj 

N0(mx)N0(x)  "  J0(x)N0(mx;J 

(4-31) 

(4-32) 


(4-33) 


from  which  we  have  an  additional  relation: 


P2-  €r,ko-h? 


^r2  ko 


e  k2  t  k2 

cr3  Ko  *  kc3 


(4-34) 
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A  simultaneous  solution  of  eqs.  (4-31)  and  (4-34)  yields  the  allowed  values 
of  h2  and  kc3.  If  a  solution  for  kc3»  that  is  h^,  really  exists,  this  solution 
is  a  surface  wave,  since  its  field  decays  exponentially  in  the  radial  direction 
and  propagates  in  the  z-di recti  on. 

By  using  graphical  methods  the  solutions  for  h^  and  k.3  are  not  so  straight 
forward.  This  is  mainly  due  to  the  fact  that  eqs.  (4-31)  can  not  be  readily 
plotted  in  the  h2b-kc3b,  (that  is  the  x-y)  plane.  Unless  such  a  plot  is  avail¬ 
able,  graphical  solution  of  eqs.  (4-31)  and  (4-34)  is  not  possible.  In  what 
follows  we  shall  proceed,  first,  to  secure  the  hgb-k  3b  plot  from  eq.  (4-31) 
by  noting  that  the  left-hand  side  of  eq.  (4-31)  is  a  function  only  of  y  (=kC3b) 
alone,  whereas  the  right-hand  side  is  function  of  both  y  and  x  (=h2b),  with  m  as 
the  parameter.  Graphical  solution  of  eq.  (4-31)  is  not  possible;  nevertheless, 
it  can  still  be  solved  by  computer. 

First,  let  us  proceed  to  simplify  eq.  (4-31)  to  a  form  that  is  more  suitable 

for  computer  solution.  This  can  be  done  by  involving  the  Bessel  function  rela- 

(251 

tionsv  '  as  follows: 

(a)  The  left-hand  side  of  eq.  (4-31)  can  be  written  with  the  aid  of  the 

relation  between  the  Hankel  function  and  the  Modified  Bessel  function 

of  the  second  kind,  K  , 

n  ’ 


Kn(z)  =  irre  i(n+1)nj  (z  ej'n)  (4-35) 

If  we  denote  the  left-hand  side  of  eq.  (4-31)  by  the  symbol  "^(y)", 
then,  by  eq.  (4-28),  we  have 

My)  ■  (jy)  Ho>>  ■  y  K°’y|  (4-36) 

where  and  K,  are  the  Modified  Bessel  function  of  the  second  kind  of 
o  1 

order  zero  and  one  respectively. 


-  43  - 


(b)  Since  the  coefficients  of  the  second  bracket  terms  both  in  the  numerator 
and  denominator  are  functions  of  x  and  y  as  well  as  the  parameters  h^  and 
h^*  we  can  simplify  this  coefficient  by,  first,  eliminating  the  parameter 
h-|,  by  the  following  consideration: 

If  we  relate  h^  and  h^  by  eq.  (4-33)  as  follows: 


hi  •  /  (  6, 


6r3*  ko  +  h3 


and  we  assume  that  region  I  is  an  empty  space  also,  then,  by  setting 
6ri  =  £r3>  we  find  hj  =  h3.  The  argument  of  the  two  Bessel  functions 

JQ  and  becomes  (jmy). 

Next,  by  eq.  (4-28)  and  by  the  relation  between  the  Bessel  function  of  the 
first  kind  with  complex  argument  and  the  Modified  Bessel  function  of  the  first 
kfnd<25),  In(z), 


I„(z)  -  e  J„(z 

we  obtain,  by  denoting  the  result  by  the  symbol  "G(y)"  as  follows: 
G(y)  =  (i)  Jo 

JQ  (Jmy^hj-l) 

-  (i)  J0(Jlny) 

J0(jmy) 

-  (- 1)  h (my> 

I0  (my) 


(4-37) 


(4-38) 


With  these  mathematical  manipulations,  we  can  write,  by  employing  the  symbol 


"fR(m,x,y)",  the  right-hand  side  of  eq.  (4-31)  as  follows: 
fn(m,x,y)  =  x  NA  +  x  G(y)  NB 


(4-39) 
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where 


Na  =  J0(X)  N^(mx)  -  J^(mx)  Nq(x) 

(4- 39a) 

Nb  =  JQ(x)  NQ(mx)  -  JQ(mx)  Nq(x) 

(4-39b) 

°A  =  Vx>  N0(mx)  '  Jo(mx)  No(x> 

(4-39c) 

db  =  J0(mx)  N0(x)  -  J0(x)  NQ(mx) 

(4-39d) 

and  G(y)  is  given  by  eq.  (4-38). 

Thus,  we  can  express  eq.  (4-31)  in  a  concise  form  as  follows: 

fL(y)  =  fR(m,x,y)  (4-40) 

where  f^(y)  is  given  by  eq.  (4-36)  and  fR(m,x,y)  is  given  by  eqs.  (4-39).  Eq.  (4-40) 
can  now  be  solved  by  computer.  The  technique  of  computer  solution  involves  the 
following  steps: 

1)  We  assign  a  suitable  value  for  the  parameter  m,  which  is  the  ratio  of 
the  inner  to  outer  radii  of  the  dielectric  tube,  hence  m  is  an  indicator 
of  the  thickness  of  the  dielectric  tube  (the  larger  the  value  of  m  the 
thinner  is  the  dielectric  tube).  To  start  with  let  m  =  0.8.  This  is 

a  case  of  thin  tube. 

2)  The  variable  y  (  *  k^b)  is  chosen  as  the  independent  variable,  which, 
for  convenience,  is  assumed  to  vary  from  0.2  to  40  in  steps  of  0.2.  We 
can  start  the  computer  by  assigning  y  =  0.2  and  f L (y )  is  computed  from 
eq.  (4-36).  This  will  give  us  a  single  positive  real  number,  say 
fL(0.2). 

3)  With  y  (=0.2)  specified  by  step  (2),  and  m  (=0.8)  chosen  in  step  (1), 
the  function  fR(0.8,x,0.2)  can  be  computed  from  eqs.  (4-39)  with  x  as 
a  dependent  running  variable,  which,  again  for  convenience,  may  be 
assumed  to  vary  from  0.2  to  40  in  steps  of  0.2.  The  result  is  tabulated, 
and  from  which  we  can  find,  graphically,  a  set  of  discrete  values  of  the 
x's  that  matches  the  following  equality  condition: 
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f L ( 0  -  2 )  =  fR(O.8,x.,0.2) 

Hence,  the  partial  solution  of  eq.  (4-31)  for  m  =  0.8  and  y  =  0.2  is  x  =  x^,  where 
x..  is  the  discrete  values  of  x  that  match  the  above  equality  condition.  It  is  to 
be  noted  that  since  both  the  modified  Bessel  functions  of  the  second  kind,  KQ(y) 
and  K^(y),  are  positive  real  numbers  for  a  positive  real  value  of  y,  hence  fL (y ) 
is  always  a  positive  real  number. 

4)  Steps  2  and  3  are  repeated  for  a  new  assigned  value  of  y,  say  y  =  0.4. 

The  matching  process  of  step  3  will  give  us  a  second  set  of  discrete 
values  of  the  x's,  say  x.,  that  will  match  the  following  equality  con- 

J 

dition: 

f L ( 0 . 4 )  =  fR(0.8,x,0.4) 

This  will  provide  us  with  a  partial  solution  of  eq.  (4-31)  for  m  =  0.8 
and  y  =  0.4,  as  x  =  x..  Obviously,  the  set  of  x.  values  will  not  be 

vJ  J 

the  same  as  x^  obtained  for  the  case  of  y  =  0.2  This  step  can  be 
repeated  and  repeated  in  steps  of  0.2  for  y,  until  y  =  40  is  reached. 

This  final  matching  process  will  give  us  the  final  set  of  partial 

solution  of  eq.  (4-31)  for  m  =  0.8  and  y  =  40,  as  x  =  x^  where  x^ 

is  a  set  of  discrete  values  of  x  that  match  the  equality  condition: 
f  L ( 40 )  =  fR(0.8,x,40). 

This  completes  the  solution  of  eq.  (4-31)  for  m  =  0.8,  with  both  x  and 
y  restrict  to  the  range  of  0  x  (or  y)C40.  The  complete  result  is  plotted  in 

the  x-y,  (that  is,  bgb-k^b),  P^ane»  shown  in  fig.  11.  A  similar  plot  for  a 

different  m  can  be  obtained  if  these  steps  are  repeated  with  different  value  of 
m.  Figs.  12  and  13  are  these  plots  for  m  =  0.5  and  0.2  respectively. 

Let  us  examine  these  graphs  more  closely  by  first  noticing  that  the  function 
fR(m,x,y),  which  is  given  in  a  precise  form  by  eqs.  (4-39),  is  discontinuous  each 
time  *  x  G(y)  Dg,  and  so  an  infinity  of  branches  occur  for  ever  increasing  x. 


¥ 
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It  passes  through  zero  each  time  the  numerator  equal  to  zero,  that  is,  =  -x  G(y)N 
The  function  f^(y),  which  is  given  by  eq.  (4-36),  approaches  zero  as  y— >  0,  and 
increasing  positively  for  increasing  y.  It  is  evident  that  for  every  y  an  infinity 
of  values  of  x  exist  which  satisfy  eq.  (4-31).  In  order  to  obtain  a  finite  number 
of  unique  solutions  of  the  mode  equation,  a  second  relation  between  x  and  y,  that 
is  between  h2  and  kcg,  is  needed.  This  is  provided,  as  mentioned  earlier  from  the 
definition  equations  (4-33)  or  eq.  (4-34), 


h2  -kc3 


*  (  6, 


r2 


€  )  Vr 
>3;  Ko 


(4-34) 


where 


kQ  =  CJ(  (UQ  £0)i  =  cj/c  =  2  7T  /  *Q,  is  the  free  space  propagation 


factor  and  aq  is  the  free  space  wavelength  correspondent  to  the  excitation 

? 

frequency  cj  .  Multiplying  by  b  yields 
1 2  .  u\2 


or 


(h2b)  +  (kc3b)  *  (  er2  ‘  £,3)  (kQb) 


x2  +  y2  =  R2 


(4-41) 


where 

R  =  (  *r2  "  *r3}*  (kob)  =  )  JTr 2  “  er3 

^0 

Eq.  (4-41)  provides  us  the  second  relation  needed  to  obtain  an  unique  solution  of 
the  mode  equation. 

We  observe  that  eq.  (4-41)  is  the  equation  of  circle  of  radius  R  centered 
at  the  origin.  The  radius  is  a  function  inversely  proportional  to  the  free 
space  wavelength  Aq.  The  graphical  solution  of  eqs.  (4-31)  and  (4-41)  is 
illustrated  in  fig.  14,  which  includes  one  of  the  x-y  plots  obtained  from  the 
computer  solution  of  eq.  (4-31),  and  superimposed  on  this  plot  are  circles 
representing  eq.  (4-41).  The  co-ordinates  of  the  common  points  of  intersection 
yield  the  values  of  x  and  y  which  fulfill  both  eqs.  (4-31)  and  (4-41).  These 
values  of  the  x's  are  the  allowed  solutions  or  eigenvalues,  h?1  =  Xj/b,  h22=  x2/b. 


20  Squares  fo  the  inch 


.  h0„  =  x  /b.  In  this  graph  the  dielectric  tube  is  assumed  to  be  made 

2p  P 

of  polystyrene,  with  €r  =  2.56,  and  the  ratio  a/b  is  taken  as  0.8.  R  is  chosen 
arbitrarily  to  be  22,  15,  and  3. 3, corresponding  to  the  ratio  of  b/  ao  =  2.8, 

1.4,  and  0.4?  respectively. 

Consider  first  for  R  =  22.  The  circle  intersects  both  branches  of  the  mode, 
or  characteristic  curves  at  P( 1 0 . 7 ,  19.3)  on  the  first  branch  curve  and  Q(20,  9.3) 
on  the  second.  The  co-ordinates  of  these  points  of  intersection  satisfy  both 
eqs.  (4-31)  and  (4-41),  and  the  x  co-ordinate  value  of  P  and  Q  gives  the  allowed 
solutions  or  eigenvalues  corresponding  to  the  ratio  of  b/  A0  =  2.8,  for  a  di¬ 
electric  tube  waveguide  having  a  dielectric  constant  €r  =  2.56.  It  is  to  be 
noted  that,  under  such  a  condition,  only  two  solutions  are  possible,  corresponding 
to  the  two  points  of  intersections  at  P  and  Q.  The  former  corresponds  to  the 
first  branch  of  the  mode  curve,  the  latter  the  second  branch.  These  are  corres¬ 
pondent  to  the  circularly-symmetric  TEg-j  and  TEq2  modes.  The  second  subscript 
"p"  on  TEgp  indicates  the  mode  number  arranged  in  order  of  increasing  cut-off 
frequency  which  is  to  be  defined  immediately.  This  is  one  of  the  feature  of 
surface-wave  guides,  that  is,  at  any  particular  frequency,  only  a  finite  number 
modes  exists.  _ 

Next,  R  =  3,3,  since  from  eq.  (4-41)  x  =  J -  y^  R.  Fig.  14  shows 
that  eq.  (4-31)  can  not  be  satisfied  for  x  3.3, that  is,  no  solution  of  the 
mode  equation  exists  for  R<  3.3.  Physically,  this  means  that  the  dielectric 
tube  waveguide  can  not  propagate  a  surface  wave  of  the  TEg1  mode,  if  the  ratio 
b/  <  0.4?.  In  other  words,  for  given  tube  dimensions  and  material,  if  aq 
is  greater  than  b/0.42,no  propagation  of  the  circularly-symmetric  TEQ1  wave 
takes  place. 

Finally,  consider  the  case, 3.3  <  R  ^  16.3a  single  unique  solution  (x,y) 
of  eqs.  (4-31)  and  (4-41)  results.  This  is  shown  in  fig.  14.  Thus,  when  R  is 
restricted  to  the  range  3.3  <  R  <  16.3  the  dielectric  tube  waveguide  which 


— 


-  52  - 


has  €r  =  2.56  and  a/b  =  0.8,  propagates  a  single  surface  wave  which  is  the 
lowest  ci rcularly-symmetric  transverse  electric  mode,  known  as  the  TEQ1  mode 
propagating  as  exp.  (j(c*>t  -  pz)). 

It  is  to  be  noted  that  as  the  dielectric  tube  wall  becomes  thicker  (i.e., 
smaller  m),  single  mode  operation  will  be  limited  to  the  longer  wavelengths. 
This  is  evident  by  inspection  of  figs.  (4-12)  and  (4-13). 

With  pairs  of  values  of  x  and  y  determined,  the  wavelength  of  the  guided 
wave  may  be  obtained  from  eqs.  (4-41)  and  (4-33). 

From  eq.  (4-41)  we  find 


f  2  2 

xfc  +  y 


(4-42) 


f  =  _c_ 


'  2  2 
x  +  y 


€r2  "  ^r3 


(4-43) 


where  c  is  the  velocity  of  light. 

Next,  from  eq.  (4-33)  we  find 

fi  =  (  UL  )  -  _JL  (  ^3 

>g  b  \  er2"^r3 

from  which  we  obtain  with  the  aid  of  eq.  (4-42) 


<»/2*  >; 

(w  a0v 


x2  + ;  er2  '  / 


(4-44) 


where  is  the  guide  wavelength. 

Eq.  (4-44)  is  shown  graphically  in  fig.  15  for  three  values  of  m.  The  cut¬ 
off  behavior  is  evident  for  all  three  cases.  Under  cut-off  condition,  for  the 
first  TEnl  mode,  with  y  »  0,  we  have  from  eq.  (4-43) 


(4-45) 


2  n  b 


r 

er2  '  €r3 


20  Sq 


The  frequency  given  by  eq.  (4-45)  is  the  lowest  frequency,  known  as  the  critical 
frequency,  for  which  a  plane  propagating  wave  may  exist  in  a  dielectric  waveguide, 
and  under  this  condition  the  wave  travels  along  the  waveguide  with  a  velocity 
characteristic  of  medium  3,  that  is,  the  medium  external  to  the  waveguide. 
Physically,  this  means  that  as  y  (  =  kc3b)  approaches  zero,  H^(jkc3b)  does 
not  attenuate  rapidly  enough  with  r  until  r  approaches  infinity,  and,  hence, 
most  of  the  transmitted  energy  through  the  dielectric  tube  waveguide  is  leaked 
out  to  medium  2.  The  characteristics  of  the  wave  propagation  along  the  wave¬ 
guide  under  this  condition  should  therefore  be  determined  virtually  by  medium 
3. 

Hence,  for  the  axial  circularly-symmetric  transverse  electric  propagating 
mode  along  the  dielectric  tube  waveguide,  there  is  a  limiting  value  of  b/  >0. 

Such  mode  of  propagation  along  the  tube  is  impossible  without  radial  dissipation 
of  energy  (radiation)  for  the  tube  with  an  outer  radius  less  than  a  limiting 
radius,  and  for  a  wavelength  longer  than  a  corresponding  limiting  wavelength 
depending  on  the  dielectric  constant  and  the  inner  radius  of  the  tube.  Immedi¬ 
ately  above  this  limiting  value  of  b /  the  wavelength  of  the  guided  wave  is 

almost  equal  to  the  wavelengths  in  vacuum  of  the  exciting  oscillation, 

With  greater  values  of  b/  the  wavelength  of  the  guided  wave  becomes  smaller 
and  approaches  the  values  ^ 

As  the  argument  of  the  Hankel  function  increases,  H^(jyr/b),  decreases 
rapidly,  and  the  energy  transfer  tends  rapidly  to  become  totally  within  the 
dielectric  tube.  Power  transfer  calculation  to  be  given  in  the  next  section 
will  validify  this.  In  order  for  the  field  of  the  propagating  waves  to  be  con¬ 
fined  to  a  thin  film  just  ouside  the  dielectric  tube  waveguide  and  stay  largely 
inside  the  dielectric  tube,  kc3  must  be  chosen  as  large  as  possible.  In  that 
case  very  little  energy  is  leaked  out  in  medium  3. 


w 
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4d.  Calculation  of  Power  Flow 

The  value  x  and  y  for  which  eqs.  (4-29)  to  (4-30)  have  physically  mean¬ 
ingful  solutions  having  been  determined  in  the  last  section,  the  coefficients, 

A,  B,  C,  and  F  may  now  be  found.  Since  for  values  of  x  and  y  which  satisfy 
both  eqs.  (4-31)  and  (4-41)  are  no  longer  linearly  independent  equations,  and 
hence  only  the  ratios  of  the  coefficients  can  be  obtained  from  eqs.  (4-29)  to 
(4-30)  directly.  Thus,  we  may  arbitrarily  choose  to  relate  the  constants,  B, 

C,  and  F  in  terms  of  the  constant  A,  then,  we  can  solve  eqs.  (4-29)  and  (4-30), 
and  eqs.  (4-29)  and  (4-30 )in  pairs  to  yield 


B 

H^(jy)  J^x)  +  j(x/y)  )( jy)  JQ(x) 

A 

H^(jy)  N-| (x)  +  j(x/y)  Hj1}(jy)  NQ(x) 

=  J0(my)  J,(x)  +  j (x/y)  Jj (my)  JQ(x) 

J0(my)  N-j  (x)  +  j  (x/y)  J^my)  NQ(x) 

(4-46a) 

C 

A 

=  JQ(x)  N-j  (x)  -  (x)  Nq(x) 

(4-46b) 

H^UyjN^x)  +  j(x/y)  H1(1)(jy)N0(x) 

F 

A 

*  J0(x)  N-|  (x)  -  J^x)  NQ(x) 

(4-46c) 

O0(my)  N-| (x)  +  J(x/y)  J^my)  NQ(x) 

where  x  and  y  are  solutions  which  satisfy  eqs.  (4-31)  and  (4-41).  As  soon  as 
the  parameters,  m,  6r2*  and  ^3  are  specified  and  the  wavelength  given, 
they  are  readily  obtained  from  the  graphical  solution  discussed  in  the  last 
section.  Hence,  the  numerical  values  of  the  ratios  B/A,  C/A,  and  F/A  can  be 
computed  from  eqs.  (46),  and  the  field  equations,  (4-24)  to  (4-26)  can  now  be 
expressed  in  terms  of  a  single  constant  which  was  arbitrarily  chosen  to  be  A. 
This  constant  will  be  determined  by  the  system  and  the  exciting  condition. 
Calculation  of  power  flow  in  the  waveguide  system  can  now  proceed  as  follows: 
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In  a  waveguide  system  employing  cylindrical  structure  of  circular  cross- 
section,  the  wave  is  travelled  along  the  z-di recti  on.  The  power  flow  is  cal¬ 


culated  from  the  general  formula: 


'i1 


l  s" ds 


(4-47) 


where 


Sn  =  J  Re  (E  x  H  )  '  n 


(4-47a) 


The  symbol  "Re"  means  the  Real  part  of  the  quantity  that  follows,  the  asterisk 
implies  the  complex  conjugate,  n  is  the  unit  vector  the  direction  of  which 
specifies  the  direction  of  the  power  flow,  and  dS  is  the  element  of  the  area 
through  which  the  power  is  flowing  out.  Total  power  flux  flowing  out  is 


calculated  by  taking  the  closed  surface  integral 


"  fs  ■ 


of  the  normal  component 


of  the  Poynting's  vector  Sn  over  the  surface  of  interest. 

The  use  of  eqs.  (4-47)  can  be  best  illustrated  by  first  considering  the 
power  flow  in  the  direction  perpendicular  to  the  outer  surface  of  the  dielectric 
tube.  The  component,  Sr,  of  the  Poynting  vector  normal  to  the  cylindrical  sur¬ 


face  of  the  dielectric  tube  can  be  calculated  as  follows: 


Sr  =  S3  *  \  =  i  Re  (  E3  x  H*)  '  ar  =  i  Re  E 


★ 

d3  Hz3 


(4-48) 


where  and  are  given  by  eqs.  (4-26),  and  when  these  equations  are  sub¬ 
stituted,  we  obtain  from  eq.  (4-47a) 

S3r.lReCC*(^)  [«(')  <h3r)]  [h*1*  Oy)]  *  (4-49) 

Using  the  following  Bessel  function  relations 


H<])  (z)  =  -  h{1}  (z) 

(«<”  ‘z>)*  ■  |>>  *  j  vz>]  * 


(4-50a) 


Jn(z)*  -  j  N  (z*) 


o'  ' 

=  H<2)  (z*) 


(4-50b) 
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.(*> 


(2 


-in  v 

eJ  )  =  - 


JnH 


«(’>  (z) 


H^Ux)  -  -  (  )  K,(x) 

H*'1  (jx)  -  -  j  (  ^-  )  K0(x) 
we  obtain  from  eq.  (4-49) 


(4-50c) 

(4-50d) 


S3r  ■  M  )2  Re  j  |c|;  (  ^  )  K0(kc3r)  K, (kc3r)  (4-51) 

=  0 

Thus,  the  component  of  the  Poynting  vector  normal  to  the  cylindrical  surface, 
S.jr»  is  pure  imaginary  and  so  there  is  no  energy  flow  out  of  the  dielectric 
tube  in  the  radial  direction.  Therefore,  the  modes  that  have  been  described 
can  be  called  surface  waves,  i.e.,  waves  which  propagate  along  the  interface 
been  two  different  mediums  without  any  transfer  of  energy  across  the  dividing 
surface. 

Next,  let  us  consider  separately  the  power  flow  in  the  direction  of  the 
wave  in  three  regions: 

(a)  Region  I,  0^  r  <a. 

In  this  region,  which  is  the  interior  to  the  dielectric  tube  where 
the  plasma  tube  is  housed,  the  Poynting  vector,  S ^ ,  is  given  by 
St  =  i  Re  E1  x  H* 

•  4  Re  E01  H*.|  (  -az)  (4-52) 

where  and  are  given  by  eqs.  (4-24).  Substituting  these 
equations  into  eq.  (4-52),  we  obtain  from  eqs.  (4-47) 

PZ1  *  i  Re  j  Sj  *  az  r  dr  d«l 

■nReF|:*  <  I  r  [vv;]  [Jo<V>]  dr 


(4-53) 
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The  two  Bessel  functions  in  the  integrand  can  be  transformed  into 

2 

(4-53a) 


M  [VV>]  =  [Il(kc3r)] 


using  the  following  Bessel  function  relations 


0o(z) 


l„U) 


-  J,(z) 


e  ■’-J*"’'  J  (z  ) 


(4-54a) 


(4-54b) 


where  In ( z )  is  the  modified  Bessel  function  of  the  first  kind  of 
order  n  with  argument  z. 

Substituting  eq.  (4-53a)  into  eq.  (4-53)  and  with  the  aid  of  the 
following  Bessel  function  relations 

j  z  I2(kz)  dz  =  -  i  z2  [l>)  -  (1  +  ^  )  I„(k z)  (4- 54c) 


and 

z  In(z)  *  -  n  In(z)  +  z  In-1(z) 
we  obtain 


=  ITa_  ( 
2 


oP 


j—)  IJ(my) 

(c3 


1 


(4-54d) 


^(my)  i((my) 


Re  FF 


(4-55) 


where  the  relation  h^ 


is  assumed. 


(b)  Region  II,  a  <  r  <  b. 

This  is  the  space  inside  the  dielectric  tube  where  the  field  is  to 
be  confined,  and  the  propagating  waves  travel  along  in  the  z-di recti  on. 
The  Poynting  vector,  Sg,  is  calculated  using  the  field  equations  given 
by  eqs.  (4-25).  The  result  is 


1 


J 


12  =  ^  Re  -2  X  -2 


■  *  Re  <-®2>  e,(2  Hr2 

=  2  i  Re  {  J^o  3 
Z  o 


)j  AA  J0(h2r)  Jq  (h2r) 

+  AB*  Jo(h2r)  N**(h2r) 

+  A*B  NQ*(h2r)  J0(h2r) 

+  BB*  NQ(h2r)  NQ*(h2r) 


(4-56) 


Substituting  eq.  (4-56)  into  eq.  (4-47),  and  using  the  Bessel  function 

I 

relation  ZQ(z)  =  -  Z-j(z),  where  the  Z's  are  cylindrical  functions,  we  obtain 

IT  ,  6^  ft  *  *  *  *  *  ★  *  ★ 

Pz2  =  Re  TT  (  -/-.p-  )  r  [AA  J]J]  +  AB  J]N]  +  BA  +  BB  N^] 

h9  J 

L  a  (4-57) 

where  for  conciseness  the  argument,  h2r,  of  all  the  Bessel  functions  has  not 
been  included  in  eq.  (4-57).  The  four  definite  integrals  of  (4-57)  can  be 
evaluated  with  the  aid  of  the  fol lowing  Bessel  function  integral  formulas 

z 

J  Zn(kz)Zn(kz)  zdz  =  iz2  ■  2Zn(kz)  Zn(kz)  -  Znl (kz)Zn+] (kz)  - 


-Zn+i(kZ)Zn.i(kz)  j 


(4- 58a) 


(the  bar  signifying  that  the  two  cylinder  functions  are  different.) 


J  [Z2(kz)]zdz  =  iz2  j  Z^2(kz)  +  (  1  -  )  Z2(kz!  |  (4- 


the  result  is  as  follows: 


AA*  J  r  [^(ty)  J*(h2r)]  dr 

=  AAV_  jj2(x)  -  J0(x)J2(x)  -  m2  [J2(mx)  -  J0(mx)J2(mx)] 


Q,  b2  A2 


(4-59a) 
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AB 


r 


r[J1(h2r)  N1(h2r)]  dr 

=  AB*b2  |  2  J1(x)N1(x)  -  Jq(x)N2(x)  -  J2(x)Nq(x) 

-  m2[2J1 (mx)N1 (mx)  -  Jo(mx)N2(mx)  -  J2(mx)No(mx)] 

2  2 
=  q2  0  * 


(4-59b) 


BA 


and 


BB 


t 


r[N1(h2r)  J](h2r)]  dr 

=  BA*b2  |  2J1(x)N](x)  -  Nq(x)J2(x)  -  N2(x)Jq(x) 

-  m2  [2J1 (mx)N1 (mx)  -  NQ(mx)J2(mx)  -  J2(mx)No(mx) ] 

"  ^3  b2  ^ 


r[N,(h2r)  N;(h2r)]  dr 
■  BB*b2 


I 


(4-59c) 


2  2 
*  Q4  *  a" 


N2(x)  -  N^(x)N2(x)  -  m2[N2(mx)  -  No(mx)N2(mx)] 

(4-59d) 


where  x  ~  h2b,  m  =  a/b,  and  ,  Q2>  Q^»  and  Q4  are  numerical  constants.  The 
purpose  of  adopting  these  numerical  constants  will  be  discussed  presently.  In 
eqs.  (4-59a)  to  (4-59d)  the  numerical  value  of  all  the  Bessel  functions,  the 

/  o  &  \ 

J's  and  N's,  can  be  obtained  either  from  tables  of  Bessel  functions'  ,  or 

from  the  subroutine  programs  of  the  computer  when  x  and  m  are  given.  With  B 

expressed  in  terms  of  A  (as  obtained  from  eq.  (4-46a)  each  of  these  definite 

integrals  eqs.  (4-59a)  to  4-59d)  will  give  a  numerical  number  in  the  form  of 
2  2 

Q  (b  A  )  as  shown,  where  Q  is  a  constant.  When  these  numerical  data  are 
substituted  into  eq.  (4-57)  we  obtain 

Pz2  =  *  (  -  )  (Qi  +  q2  +  Q3  +  Q4)  (b2  A2) 


> _ 


(4-60) 


where,  as  mentioned  above,  for  simplicity  in  expression,  ,  Q,,,  Q^,  and  are 
numerical  numbers  computed  from  eqs.(4-59a)  to  4-59d)  in  conjunction  with  eq. 
(4-46a) . 


(c) 


Region  III,  r  >  b. 

This  is  the  region  of  space  surrounding  the  dielectric  tube,  which  has 
been  assumed  to  be  the  empty  free  space  in  previous  analysis.  The 
Poynting  vector,  Sj,  is  calculated  using  the  field  equations,  eqs.  (4-26). 
The  result  is 


h 


=  i  Re  Eg  x 
=  -  az  i  Re  E^3  Hr3 

=  a  i  Re  CC*  (  -J  ^  \  -  )  [H^1 }  (h^HtH*1 } '  (h,r)  ]* 

h3  h3 


(4-61) 

The  power  flow  is  calculated  in  accordance  with  eq.  (4-47)  and  (4-47a). 
Substituting  eq.  (4-61)  into  (4-47)  yields 

*  .  4) 


z3 


>1  Re  CC  (  ) 


I 


°°  r[(H^1)’(h3r)]  [H^'Oy)]*  dr 

(4-62) 


The  integral  can  be  evaluated  by  first  using  the  Bessel  function  relations 
of  eqs.  (4-50a)  to  (4-50d)  to  transform  the  integrand 

[n©1  > ' (h3r)]  [H^'(h3r^  =  )  Kl(kc3r^ 

then,  using  the  Modified  Bessel  function  integral  formula 

j  z[K2(kz)]dz  =  iz2  {  K*2(kz)  -  (1  +  ^  )  K2(kz )  j  (4-63a) 
and  the  recurrence  formula 

z  K’(z)  =  -  n  Kn(z)  -  z  Kn-1 (z)  (4-63b) 

to  evaluate  the  integral.  When  the  result  is  substituted  into 
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eq.  (4-62)  we  obtain 

p23  ■  |c|2  (  y>  b4  (  )  K(<*> 

(4-64) 

The  magnitude  of  power  flow  carried  by  the  traveling  waves  in  the  three 
regions  can  be  best  illustrated  by  the  following  example: 

With  reference  to  fig.  14,  the  solution  that  satisfies  eqs.  (4-31)  and  (4-41) 
for  the  polystyrene  (  =  2.56)  dielectric  tube  waveguide  having  an  a/b  ratio  of 

0.8,  is 

x  =  9.28  and  y  =  11.78. 

The  waveguide  is  operating  in  the  TEg-j  mode  at  a  frequency  corresponding  to  the 
b/  =  1,9  (that  is  R  =  15). 

With  x,  y,  m,  and  b/  given,  the  value  of  mx  and  my  can  be  computed  and 
all  Bessel  functions  needed  to  compute  the  power  flow  can  be  obtained  either  from 
the  tables  of  Bessel  functions or  from  the  subroutine  programs  of  the  computer. 
We  can  proceed  first  to  compute  the  ratios,  B/A,  C/A,  and  F/A  from  eqs.  (4-46a), 

(4- 46b )  and  (4-46c)  respectively  to  yield: 

B  =  -  0.23055  A 
C  »  -  j  0.873  e11 ,78  A 
F  =  -  0.60577  e-9-424  A 

Next,  the  terms  (  c)jjq  p)  which  appears  in  all  power  flow  formulas,  can  be 
computed  from  the  data  b/  ;*  =  1.9  and  the  following  relation 

h3  ■  «kc3>2  *  *r3  ko  -  J 
given  by  eq.  (4-27c).  The  procedure  is 

- (  *V  V  >3  (  (  6  >  '  (  “6  > 

=  6.8714  x  104  b"2. 


nf(y>  y  ( 


With  these  data  available,  the  magnitude  of  the  power  flow  in  the  three 
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regions  can  be  computed  from  eqs.  (4-55),  (4-60),  and  (4-64)  to  yield: 

Pzl  =  0.15  e'18’848  (  fra2)  A2 
Pz2  =  9.62  (71b2)  A2 
Pz3  =  0.93  (71b2)  A2 

The  value  of  Pz2  is  computed  based  on  the  value  of  Q1 ,  Q2,  Q3>  and  Q4  given  by 
eqs.  (4-59a)  to  (4-59d).  They  are 
Q1  =  3.22  x  10"3 


Q2  =  0.3  x  10 
Q3  =  0.24  x  10 
=  0.06  x  10 


-3 


-3 

-3 


By  comparing  the  relative  magnitudes  of  P^ ,  Pz2>  and  Pz3,  we  conclude  that: 

(1)  Power  carried  by  the  wave  enclosed  within  the  inner  radius  of  the 
dielectric  tube  is  practically  zero. 

(2)  The  power  carried  by  the  wave  exterior  to  the  dielectric  tube  is 
about  9.6%  of  Pz2> 

Further  detailed  calculations  of  Pz3  as  function  of  distance  from  the  surface  of 
the  dielectric  tube  show  that  almost  99%  of  Pz3  are  confined  within  a  region  of 
space  having  a  radius  between  b  and  1.2b,  corresponding  to  a  distance  from  r  =  b 
to  r  =  1.2b  =  1.2  x  1.9  =  2.28  P\Q.  (see  fig.  16).  At  this  distance,  the 

relative  electric  field  strength  as  compared  with  its  value  at  r  =  b  can  be 
computed  as  follows: 

-  ^'(I.Z  jkc3b) 

Hr'Okc3b) 

.  K,  (t.2y) 
k,-  Ty) 

.  Kj (14.136) 

~  x^rra) 


'E""(r=bj” 


=  0.002 


This  implies  that  in  a  distance  of  over  two  wavelengths  the  electric  field  has 
decayed  to  a  relatively  small  value  about  0.2%  from  its  value  at  r  =  b. 

V.  The  Cylindrical  Metal  Waveguide  Partially  Filled  with  a  Coaxial  Dielectric 


Tube. 

With  reference  to  fig.  17,  a  high  dielectric  constant  dielectric  tube  of 
inner  and  outer  radii  "a"  and  "b"  respectively  is  coaxially  fitted  into  a  cylin¬ 
drical  metal  waveguide  having  an  inner  radius  b.  The  interior  of  the  metal 
waveguide  of  circular  cross-section  is,  thus,  divided  into  two  regions  of  diff¬ 
erent  dielectrics. 


Dielectric  tube 


Fig. 17.  The  cylindrical  metal  waveguide  partially 
filled  with  a  coaxial  dielectric  tab*. 


M 
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(a)  Region  I,  (0<r^a):  This  is  a  region  of  space  which  may  be  a  free 
space  or  some  other  isotropic  dielectric  medium  of  small  dielectric 
constant,  and  in  practice,  this  may  be  the  empty  space  where  the  plasma 
tube  is  to  be  housed, 

(b)  Region  II,  (a£r4b):  This  is  the  region  of  space  the  electromagnetic 
fields  and  waves  are  to  be  confined  in  and  guided.  The  dielectric 
material  used  is  assumed  to  be  isotropic  having  a  dielectric  constant, 

^r2>  higher  than  that  of  the  region  I. 

It  is  also  assumed  that  the  conductivity  of  the  two  dielectric  mediums  to 
be  zero,  and  that  of  the  metal  cylindrical  tube  to  be  infinite. 

With  these  specifications,  this  type  of  waveguide  structure  may  be  considered 
as  a  modified  form  of  the  more  general  dielectric  tube  waveguide  discussed  in  the 
last  section,  in  which  region  III  is  replaced  by  a  grounded  conducting  cylindrical 
tube  having  an  inner  radius  "b".  Furthermore,  since  we  are  also  concerned  with 
the  field  solution  of  the  circularly-syircnetric  transverse  electric  types,  we  can 
proceed  to  analyze  the  field  problem  using  the  results  already  obtained  in 
section  IV  without  going  into  further  repetition  of  the  details  in  the  solution 
of  the  boundary  value  problem.  These  results,  which  have  been  listed  in  eqs. 
(4-24)  to  (4-25),  are,  for  convenience,  reproduced  in  the  following: 

Region  I:  0^*r^a 
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H 
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Region  II:  a^r^b 


'z2  ’  [A  Jo(h2r)  *  B  W1 

(5-2a) 

o 

II 

CM 

i 

ii 

CM 

[AJo<h2r)  * 

B  N‘(h2r)] 

e'^z 

i 

(5-2b) 

r2  "  °^Po 
i 2  ■  -R“- 

A  Jo^h2r^  + 

B  N^(h2r) 

- 

e-  tz 

(5-2c) 

z2  =  0 

where 


hf  =  kf  +  V2  =  er1  k2  -  (S2  (S-3a) 

4  ■  4  *  t 2  ■  er2  k2  -  P  2  (5-3b) 

where  2f  is  identically  the  same  in  h^  and  h^.  and  for  propagation  mode,  if  58  J  f>  * 
is  the  phase  shift  constant.  K  =  6^.  =  k0*/6y* 

In  eqs.  (5-1)  to  (5-2),  A,  B,  and  F  are  new  constants  which  are  to  be  deter¬ 

mined  by  a  new  set  of  boundary  conditions  appropriate  to  this  new  type  waveguide 
structure.  The  boundary  conditions  are: 

(a)  At  the  boundary  between  two  dielectric  mediums,  both  and  H,  must  be 
continuous,  that  is, 


r=a 


EU2 

,  and  Hz^ 

Hz2 

r=a 

r=a 

(b)  At  the  wall  surface  of  the  metal  cylindrical  tube,  the  tangential  com¬ 
ponents  of  the  electric  field  must  vanish,  that  is. 


Substitution  of  these  boundary  conditions  into  eqs.  (4-1)  and  (4-2)  gives  us 
F  hj1  J^a)  *  h“J  [a  j’dy)  +  B  N^a)]  (5-4a) 

F  JQ(h1a)  =  A  J0(h2a)  +  B  NQ(h2a)  (5-4b) 
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A  J^h^)  +  B  %(h2b)  =  0  (5-4c) 

In  order  to  obtain  a  solution  different  from  0  for  A,  B,  F,  the  determinant  of 
the  coefficients  must  vanish.  Thus,  by  letting 

h2ljo(h2a)  h2lNo(h2a)  ’  hl'1jo(hla) 

Jo(h2a)  No(h2a)  ‘  Je(hla)  =  0 

Jo(h2b)  No(h2b)  0 


and  expansion  of  this  determinant  leads,  if  we  discard  trivial  solut,  to  the 
transcendental,  or  the  mode  equations 

y  J0M  =  x  J0(x)  NQ(mx)  -  JQ(mx)  NQ(x)  (5-5) 

0o(y)  J0(x)  NQ(mx)  -  JQ(mx)  NQ(x) 

where,  for  convenience  of  notation,  the  following  substitutions  have  been  made: 
m  =  b/a  >  1,  x  =  h2a,  and  y  =  h^a 

Furthermore,  using  the  recurrence  relation  of  the  Bessel  functions  Zn(z) 

Zn(z)  =-Wz)  +  (n/z)  Zn(z) 

where  Zn  is  the  cylinder  functions,  Jn(z)  or  Nn(z),  eq.  (5-S)  may  be  expressed 
in  the  form 

L  (y)  =  R  (m,x)  (5-6) 


where 


Jo(y) 

l  (y)  =  y  j-fjy 


(5-6a ) 


is  the  left-hand  side  function  of  eq.(5-5),  which  is  a  function  only  of  y  alone; 


R  (m,x) 


JQ(x)  N^mx)  -  J-|  (mx)  Nq(x) 


(5-6b) 


is  the  right-hand  side  function  of  eq.  (5-6),  which  is  a  function  of  x  with  m  as 
the  parameter. 


x  and  >  arc  further  conn 
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(5-3b),  that  is 

* 2  -  >4  -  €r2  *o  *  •>?  -  6r,  <  <^> 

hence, 

h2  '  hl  =  (  er2  ‘  6rl)  ko  (5-7a) 

where  kQ  =  (2  ft/  ^Q)  =  cO  is  the  free  space  propagation  factor  and 

o 

is  the  free  space  wavelength.  Multiplying  both  sides  of  eq.  (5-7a)  by  a 
yields 

x2-y2  =  (  €r2-  6rl)  (k0a)2 

=  (  ^2  *  £r])  (  13±  )2  (5-7b) 

A) 

which,  together  with  eqs.  (5-6),  completely  determines  the  frequency  and  parameters 
of  the  waveguide  structure. 

Graphs  of  the  functions,  L(y)  and  R(m,x),  for  three  values  of  m  are  plotted 
in  figs.  (19),  (20),  and  (21).  The  function  L(y)  is  discontinuous  each  time  J-j(y) 
vanishes  and  so  an  infinity  of  branches  occur  for  ever  increasing  y.  It  passes 
through  zero  each  time  JQ(y)  has  a  zero.  The  function  R(m,x)  vanishes  each  time 
the  denominator  of  R(m,x)  vanishes  and  so  there  are  also  infinity  of  branches 
for  ever  increasing  x.  Zeros  of  R(m,x)  occur  each  time  the  numerator  of  R(m,x) 
has  a  zero.  It  becomes  evident  from  figs.  (19)  to  (21)  that  for  every  finite  y, 
an  infinity  of  x  exist  which  satisfy  eq.  (5-6)  and  vice  versa.  The  relations 
between  x  and  y  for  three  values  of  m  are  plotted  in  figs.  (22),  (23)  and  (24). 

In  order  to  obtain  a  finite  number  of  unique  solutions  of  the  mode,  or  character¬ 
istic  equation,  (5-5)  or  (5-6),  a  second  relation  between  x  and  y  is  needed. 

This  is  provided  for  by  eq.  (5-7b),  which  is  the  equation  of  the  rectangular 
hyperbola.  The  graphical  solution  of  eq.  (5-6)  is  illustrated  in  fig.  25',  which 
includes  only  the  first  few  branches  of  the  x-y  plot  for  an  arbitrary  value  of 
m,  and  superimposed  on  them  is  the  rectangular  hyperbola  representing  eq.  (5-7b) 
for  a  given  wavelength  >Q.  We  observe  that  for  a  given  frequency,  that  is, 

for  a  given  *Q,  only  one  set  of  solutions,  called  the  mode  solutions,  (x^y.) 


. 


R  2*  70  20  -iO  Squares  to  the  inch 
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with  i  -  1 ,  2,  3,  . . .  p,  is  possible.  These  modes  are  the  TEQ1 ,  TEQ2,  ...  TE0p 
modes.  The  corresponding  wavelength  in  the  waveguide  is  found  from  the  relation 
given  by  eq.  (5-3b),  that  is 

An 


<£ 


r2 


f  x 

271a 

^  Ac 


(5-8) 


la  j 

'o  J 


where  p\  is  the  wavelength  in  the  guide.  We  note  that  as  the  denominator  of 
(5-8)  approaches  zero,  and  for  a  given  AQ»  Ag  becomes  longer  and  longer. 

In  the  limit  when 


T2 


x2 

(  ?2ta)2 
>0 


=  0 


or 


.  <Lr->  JTi 


(5-9) 


If  frequency  is  increased  (  /\Q  decreased),  x  increases  while  y  decreases  until 
at  a  critical  frequency,  f  ,  y  becomes  zero.  When  this  happens,  eq.  ( 5 -7b) 


reduces  to 

,  2  K  a, 
Xc .  ( 

which  may  be  written  as 


r2 


rl 


=  ( 


r2 


>1 


(5-10) 


where  x  and  A  are  the  value  of  x  and  A  at  the  critical  condition,  and  A 
c  c  1 

is  called  the  critical  wavelength.  Thus,  when  y  =  0,  AQ  — *>  >c,  eq.  (5-8) 
reduces  to  Ag  =  Aq/  J £r2.  Now,  as  y-»0 


J0(y) 


Lim  L(y)  =  Lim  y 

y  — ►  0  y-+  0 


►  1 


(5-11) 


then,  the  relation  of  eqs.  (5-6)  and  (5-6b)  gives  us 
JQ(x)  N^mx)  -  J-|  (mx)  Nq(x) 

1  =  X  J-j  (x)  N^mx)  -  J-j  (mx)  N^x) 


(5-12) 


k 
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Values  of  the  first  few  x  that  satisfy  eq.  (5-12)  can  be  computed  for  three  values 
of  m  are  listed  in  the  following: 

m=1.2  m=1.6  m=2.0 

2 

4.9 
7.95 
11.5 


cl 

kc2 


8.5 


3 

8 

13.2 


c4 


We  note  that  there  are  some  discrepancies  in  the  values  of  the  xc1 ,  x^... 

listed  above  as  compared  with  that  which  are  obtained  from  figs.  (5-22),  (5-23), 
and  (5-24),  mainly  due  to  the  inaccuracy  in  the  graphical  method  in  obtaining 
these  figures. 

The  results  obtained  so  far  show  relatively  little  value  to  our  objective. 
However,  if  one  attempts  to  go  into  the  matter  deeper,  one  would  find  that  a 
further  increase  in  frequencies  (that  is,  decreased  in  7*  ),  eq.  (5-7b)  can  no 
longer  be  satisfied  unless  y  becomes  imaginary.  Let  y  (  ■  h,a)  =  -  j  kc,a, 
where  (k^a)  is  the  magnitude  of  y,  a  real  positive  quantity.  The  negative  sign 
chosen  is  based  on  the  consistency  requirement  to  be  discussed  immediately. 

Now,  by  replacing  y  by  (-jk^a),  and  usinc  the  following  modified 
Bessel  function  relation 

I„(z)  -eTJJn*  Jn(z  e  ) 
the  left-hand  side  function  of  eq.  (5-6)  can  be  written  as 

.  Jol-jk-i3) 

L  (-■>kcl4>  ‘  <-Jkc,a) 


\l*> 


'Ai*> 


(5-13) 

.  .  •  r  ns  una  I  ter  i*d,  and  for  convenience  in 
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R(m,x)  =  x 


JQ(x)  N1  (mx)  -  J-,  (mx)  NQ(x) 


J 1  ( x )  N-j  (mx)  -  J-j(mx)  N-j  ( x) 

Instead  of  eq.  (5-6),  we  now  have  a  new  mode  equation 


(5-14) 


L  (kcla)  =  R  (m,x) 


(5-15) 


where  L  (k^a)  is  given  by  eq.  (5-13)  and  R(m,x)  is  given  by  eq.  (5-14). 
note  that  as  (k^a)  approaches  zero,  eq.  (5-13)  becomes 


We 


Lim  L  (kc1a)  =  Lim 


IQ  (kcla) 

<kcia> 


kcla 


kcla 


which  is  consistent  with  the  result  of  (5-11).  If  a  positive  sign  have  been  used 
for  imaginary  y,  the  limiting  value  of  (5-13)  would  have  been  (-1)  which  is  not 
in  agreement  with  eq.  (5-11). 


I 

It  is  to  be  noted  that,  unlike  L(y)  of  eq.  (5-6a),  L  (k^a)  is  always  positive 
real.  It  increases  with  increasing  (k^a),  whereas  the  function  R(m,x),  which  has 
been  discussed  in  details  earlier,  is  a  discontinuous  function  having  infinity  of 
branches  for  ever  increasing  x.  Hence,  in  order  for  eq.  (5-15)  to  be  satisfied, 
no  negative  values  of  R(m,x)  should  be  included  in  forming  the  graphical  solution 
of  eq.  (IS).  This  is  shown  in  figs.  (5-25),  (5-26),  and  (5-27).  Figs.  (5-28), 
(5-29),  and  (5-30)  are  the  x-k^a  plots  corresponding  to  figs.  (5-25),  (5-26), 
and  (5-27)  respectively. 

Finally,  in  the  limit,  when  *0,  and  ( k^-j a )  — < ►OO.  It  is  more  convenient 
to  rewrite  eq.  (5-15)  in  the  form 


L'(kcia) 


R(m,x) 


(5-13) 


-1 


Then,  by  taking  the  limit  of  [L  (kcla)]  as  (kc^a)-*oo,  we  find 


Lim 


kc,a-.oo 


kcla-*» 
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N-j  (x) 


J-|  (mx) 


N-j  (mx) 


-*■  0 


and  under  such  a  condition,  in  order  for  eq.  (5-13)  to  be  satisfied 
1  1  J^x)  N-|  (mx)  -  (mx)  N^x) 

R(m,x)  x  JQ(x)  N-j  (mx)  -  J^mx)  Nq(x) 

hence,  Jj(x)  N-j(mx)  -  J^mx)  Nj(x)  =  0.  This  gives  us 

^(x) 


(5-14) 


x  is  the  root  of  eq.  (5-14),  and  the  first  few  roots  for  three  values  of  m  are 
listed  in  the  following: 


m 

x  (  7\0— **0  ) 

1.2 

15.72 

31.43 

1.6 

5.22 

10.52 

15.73 

20.95 

31.42 

36.68 

2.0 

3.19 

6.31 

12.58 

15.72 

18.86 

22.0 

25.15 

To  complete  the  graphical  solution,  the  connecting  equation,  (5-7b)  should 
also  be  modified.  This  is  accomplished  by  substituting  (-jkcla)  for  y  in  eq. 
(5-7b)  to  yield 


X2  *  <kc,a)2  -  <  6r2  -  «r,  )  (  ^  >2 

=  R2  ,  (5-13) 

where 

*  ■  AT7" ^  >  (5-,3a) 

Eq.  (5-13)  represents  a  circle  of  radius  R  centered  at  the  origin.  Technique  of 
graphical  solution  of  eq.  (5-15)  and  (5-13)  is  illustrated  by  fig.  (5-31),  in 
which  the  (x-kcl a)  plot  is  that  which  reproduced  from  fig.  (5-28)  whereas  the 
radius  of  the  circle  of  eq.  ( 5-1 3a )  is  chosen  arbitrarily  to  be  20  and  6^ 
is  taken  to  be  2.56  and  €  j  =1. 

By  comparing  fig.  (5-31)  with  fig.  (14),  we  notice  the  similarity  between 
these  two  plots  is  remarkable,  and  Indeed,  by  a  further  study  of  the  general  feature 


of  fig.  (5-31)  to  be  described  and  discussed  presently,  one  would  find  several 
properties  and  characteristics  that  are  common  to  both  this  type  of  waveguide 
and  the  dielectric  tube  waveguide  discussed  in  section  IV. 

First,  we  note  that,  as  in  the  case  of  a  dielectric  tube  waveguide,  cut-off 
occurs  when  R  (that  is  the  ratio  of  a/  7\Q)  less  than  a  certain  critical  value  of 
x  for  which  (k^a)  =  0-  The  critical  values  of  x  are  xcl  =  8.8  for  the  TEQ1 
mode,  and  x^  =  23.7  for  the  TEq2  mode.  When  R  is  chosen  to  lie  between  xci^R<xC2« 
the  circle  intersects  only  the  first  branch  of  the  mode  characteristic  curve  at 
one  point  whose  co-ordinate  is  shown  to  be 
x  =  12.4  and  k^a  =  15.7 

for  R  =  20  as  shown  in  the  figure.  Thus,  a  single  mode  operation  is  ensured  by 
a  proper  choice  of  parameters.  When  R'>XC2>  the  circle  will  intersect  the  mode 
characteristic  curves  (two  of  such  curves  are  shown  in  the  figure)  at  more  than 
one  point  depending  upon  the  value  of  R  chosen  and  the  parameter  m  used.  Hence, 
multiple  modes  always  exist  just  like  many  other  metal  waveguides  when  it  is 
operated  at  very  short  wavelengths.  Nevertheless,  there  are  several  parameters 
that  are  readily  available  to  control  or  restrict  the  number  of  modes,  the 
existence  of  which  might  hinder  the  proper  operation  of  the  device.  For  example, 
the  dielectric  tube  may  be  made  thinner  (that  is,  smaller  m)  so  that  the  cut-off 
point,  x^,  may  be  made  to  move  upward  in  value.  This  will  extend  the  maximum 
range  of  R  (that  is  the  a/  ratio)  that  can  be  achieved  for  a  single  mode 
operation.  The  same  result  could  be  accomplished  by  a  reduction  in  the  plasma 
tube  radius  (that  is  the  radius  a  also)  if  this  is  feasible.  If  on  the  other 
hand,  an  over-sized  dielectric  tube  must  be  used,  one  still  can  attain  high 
frequencies  by  limiting  the  number  of  higher  order  modes  that  could  exist  to 
a  smaller  number. 

Fig.  5.32  is  the  (  *Q/  *g)  vs.  (a/  aq)  plot  for  fr2  =  2.56  and  6^  *  1. 

This  is  obtained  by  relating  eq.  (5-1 3a )  to  eq.  (5-3b)  using  a/*Q  as  the 
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independent  variable.  When  this  is  compared  with  the  thin  tube  curve  (m  =  0.8) 
of  the  fig. 15  the  similarity  in  the  basic  property  between  the  dielectric  tube 
waveguide  and  the  metal  waveguide  which  is  partially  filled  with  coaxial  dielectric 
tube  operating  in  the  very  short  wavelength  modes  is  further  evident. 

In  the  next  step,  we  shall  proceed  to  determine  the  three  arbitrary  constants, 
A,  B,  and  F  that  appear  in  the  field  equations  (5-1)  and  (5-2).  This  is  easily 
obtained  by  solving  eqs.  (5-4a),  (5-4b),  and  (5-4c)  simul taneously  to  yield 
B  and  F  in  terms  of  A  as  follows: 


JQ(y)  ^U)  -  (x/y)  J^y)  JQ(x) 
J0(y)  N-j  (x)  -  (x/y)  J-j  (y)  NQ(x) 
J0(x)  N-j  (x)  -  J^x)  NQ(x) 

0o(y)  N-j  ( x )  -  (x/y)  J-|  (y)  NQ(x) 


(5-15) 

(5-16) 


Consider  only  the  case,  y  =  -jkcla,  which  is  of  great  importance  to  us. 
By  substituting  f-jkcla)  for  y  and  using  the  Bessel  function  relations 

1„U)  =  eTJin"  0„U  <^‘1 

we  obtain  from  eqs.  (5-15)  and  (5-16) 

Vkcla>  J1(x)  ■  (x/kcla)  Mkcla*  Jo(k)  ,, 

B  =  -  A - (5'  7) 

I0(kcla)  Ni(x)  -  (x/kcla)  I -j ( kcl a )  NQ(x) 

3nd  J0(x)  N  (x)  -  J,(x)  N  (x) 

r  A  °°  10 _  (5-18) 

Vkcla)  Vx)  ■  (x/kcla^  X1  ^kcla^  Vx) 


where  x  and  kc^a  are  the  solution  of  the  mode  equation,  the  various  Bessel  functions 

/  or  \ 

can  be  readily  computed  either  from  the  tables  of  Bessel  function'  '  or  from 


the  subroutines  of  the  computer  program,  hence,  the  constants  B  and  F  can  be 
computed  in  terms  of  a  single  constant  A  which  depends  on  the  system  design  and 
the  experimental  conditions. 
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With  the  ratios,  B/A  and  F/A,  determined  the  electric  field  in  the  two 


regions  can  be  computed  from  eqs.  (5-lc)  and  (5-2c)  as  follows: 


•01 


■02 


-  JF  ( 


kcla 


>  Vkcla  5  > 

for  o4r<a 


v-:>  {■•!*$) 


for  a  <  r  <  b 

The  propagation  factor  e^  ^ z ^  is  assumed. 


(5-19) 


(5-20) 


Fig.  (5-33)  shows  the  electric  field  distribution  in  the  two  regions  as 
computed  from  eqs.  (5-19)  and  (5-20).  The  parameters  and  data  used  to  compute 
eqs.  (5-19)  and  (5-20)  were  taken  from  the  graphical  solution  of  the  mode  equation 
as  shown  in  fig.  5-31,  with  R  =20,  =  2.56,  =  1 ,  m  (  =  b/a)  =  1.2, 

and  x  =  12.4  and  k^a  =  15.7.  The  capability  of  this  waveguide  structure  to 
confine  the  electric  field  in  the  high  dielectric  region  is  evident. 

We  conclude  that  a  high  dielectric  constant  dielectric  tube,  when  it  is 
coaxially  tight-fitted  into  a  cylindrical  metal  tube,  represents  a  new  kind  of 
waveguide,  limited  on  the  outside  by  the  conducting  surface  of  the  metal  tube  and 
on  the  inside  by  a  medium  of  smaller  dielectric  constant.  For  high  frequencies 
the  waves  are  confined  to  the  space  occupied  by  the  high  dielectric  constant 
medium,  andjif  the  dielectric  tube  is  thin,  propagation  of  a  single  mode  at  very 
short  wavelength  becomes  a  reality. 


VI.  The  Hollow  Dielectric  Tube  Waveguide  Supported  by  a  Conducting  Cylinder. 

A  thin  dielectric  tube  supported  Internally  by  a  grounded  metal  conducting 
tube  forms  the  third  class  of  waveguide  structure  to  be  Investigated  in  this 
studies.  As  shown  in  fig.  34,  the  Inner  and  outer  radii  of  the  dielectric 

— _ . _ _ _ _ _ _ _ _ J 
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Region  III 


Grounded  conducting  cylinder 


Region  II 


Pig. 34.  The  hollow  dielectric  tube  waveguide  supported 
by  a  conducting  cylinder 

tube  sleeving  is  a  and  b  respectively.  The  outer  radius  of  the  metal  tube  is  a, 
which  is  tightly  fitted  into  the  dielectric  tube.  This  class  of  waveguides,  like 
the  dielectric  tube  waveguide  investigated  in  section  IV,  generally  called  open¬ 
boundary  waveguides,  is  characterized  by  having  the  field  outside  the  boundary. 
These  guides  support  modes  of  propagation  called  the  surface  waves,  waves  which 
propagate  along  the  interface  between  two  different  mediums  without  any  transfer 
of  energy  across  the  dividing  surface.  The  field  of  the  surface  waves  decays 
exponentially  in  the  direction  normal  to  and  away  from  the  guiding  structure, 
and  by  an  appropriate  choice  of  parameters  the  waves  and  fields  can  be  confined 
within  a  narrow  space  surrounding  the  dielectric  tube. 

The  problem  of  surface  wave  propagation  along  a  dielectric  coated  solid 
wire  and  the  practical  usefulness  of  such  a  line  for  microwave  transmission 

have  been  investigated  by  Goubau^2^,  Barlow  and  Cullen^2®),  Barlow  and  Kar- 

(291  (301  (311 

bowiak'  ,  Roberts'  .  A  more  extensive  account  may  be  found  elsewhere' 

This  type  of  waveguide  structure  possesses  some  very  desirable  properties 
as  transmission  lines  at  millimeter  and  submillimeter  wavelengths.  However, 
in  view  of  its  structural  limitations,  this  type  of  structure  is  not  suitable 
to  be  Integrated  directly  with  the  plasma  tube,  unless  a  length  of  the  support¬ 
ing  metal  tube  Is  appropriately  cut-away  to  make  room  for  the  plasma  tube.  In 
that  case  the  part  of  the  surface  waveguide,  with  its  Internal  supporting  metal 


tube  removed,  operates  as  the  dielectric  tub  waveguide  modes  discussed  in  section 
IV,  whereas  the  remaining  part  of  the  surface  waveguide  serves  purely  as  a 
microwave  transmission  line. 

The  method  of  mathematical  analysis  of  the  surface  waveguide  is  very  similar 
to  the  other  two  types  of  waveguides  already  analyzed  and  discussed  in  details 
in  sections  IV  and  V.  In  what  follows  only  a  brief  outline  of  the  mathematical 
analysis  will  be  presented  with  the  result  of  the  graphical  solution  given. 

The  presence  of  the  grounded  conducting  tube  interior  to  the  dielectric 
tube  sleeving  excludes  the  region  I  from  any  consideration  as  far  as  the  field 
solution  is  concerned.  Instead  of  solving  a  boundary  value  problem  involving 
three  boundaries  we  need  to  match  the  fields  at  two  boundaries  only.  This 
will  result  in  a  much  simplification  of  the  field  problem.  Thus,  let  us 
consider 

(a)  Region  II,  a<.r<b:  This  is  the  region  of  space  filled  with  high 
dielectric  constant  material  where  the  fields  and  waves  are  to  be 
confined  and  guided.  The  dielectric  material  used  is  assumed  to  be 
Isotropic  and  lossless. 

(b)  Region  III,  r>b:  This  is  the  region  of  space  surrounding  the  di¬ 
electric  tube  sleeving.  It  is  assumed  to  be  the  empty  free  space. 

With  these  specifications,  this  waveguide  structure  may  again  be  considered 
as  a  modified  form  of  the  more  general  dielectric  tube  waveguide  discussed  in 
section  IV,  in  which  region  I  is  replaced  by  a  grounded  conducting  cylindrical 
metal  tube  having  an  outer  radius  a.  Again,  since  we  are  concerned  only  with 
the  field  solution  of  the  circularly-synmetrlc  transverse  electric  types, 
we  can  proceed  to  analyze  the  field  problem  using  the  results  obtained  in 
section  IV  without  going  Into  further  repetition  of  details  in  the  solution 
of  the  boundary  value  problem.  These  results,  which  have  been  listed  In  eqs. 
(4-25)  to  (4-26),  are  reproduced  in  the  following: 
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Region  II :  a  4  r  <  b 

Hz2  '  [A  Jo(h2r>  *  B  N0(h2r>  ] 


H«-° 

Hr2  *  *4- 
h2 

J 


J'(h2r)  +  B  H^(h2r)'j 


Er2-° 


E«*--(i2"  [A  W>  tBNi<h2r>] 


-  #z 


Ez2  =  0 


(6-la) 


(6-lb) 


(6-lc) 


Region  III:  r  >  b 

HZ3  =  CHo',<h3r>  e'^ 

H«3  =  ° 

"r3--Ch;  e-yz 

Er3  =  ° 


Ez3  *  0 


where 


h|  ■  kf  +  I2 

h3  ■  Ukc3>2 

*  kj  +  f 


fr2ko-  I5' 


6r3k0-  P2 


(6-2a) 


(6-2b) 


(6-2c) 


(6-3a) 


(6-3b) 


/. 


where  the  h's,  k's,  ^  ,  |&  ,  have  been  defined  before,  and  the  choice  of 
(jkcj)  for  hj  for  region  III  has  also  been  discussed  in  detail  in  section  IV, 

No  further  explanation  Is  needed. 

The  constants,  A,  B,  and  C  can  be  determined  by  requiring  both  the  electric 
and  magnetic  fields  to  satisfy  the  following  boundary  conditions: 

(a)  On  the  surface  of  the  metal  cylindrical  tube,  the  tangential  components 
of  the  electric  field  must  vanish,  that  is 
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■02' 


r=a 


(b)  At  the  boundary  between  two  dielectric  mediums,  both  E^  and  Hz  must 
be  continuous,  that  is 


l02 


=  E 


03 


,  and  HzJ 


=  H 


z3 


lr*b 


lr»b 


lr*b 


r  ab 


and  from  which  we  obtain 

A  0o(h2a)  +  B  N'(h2a)  =  0 


h*1  A  0'(h2b)  +  B  N'(h2b)  -  h’1 


(1) 


A  JQ(h2b)  +  B  NQ(h2b)  -  C  (h3b)  =  0 


1 

C  H*1 }  (h3b)  =  0 

(6-4a) 

(6-4b) 

b)  =  0 

(6-4c) 

In  order  to  obtain  a  solution  different  from  zero  for  A,  B,  and  C,  the  deter¬ 
minant  of  the  coefficients  must  vanish.  Thus,  by  letting 


J^(h2a)  N;(h2a)  0 

i.-'  j'  (h2b>  h’1  N'(h2b)  -  h;1  h<'>'  (h3b> 


•  H<')(h3b) 


=  0 


J0(h2b)  N0(h2b) 
and  expansion  of  this  determinant  leads  to  the  transcendental,  or  mode  equations 

(6-5) 


(Jy)  H°'’,  (Jy) 

TO1 


H 


(jy) 


O0(mx)  N  (x)  -  J  (x)  N'  (mx) 
X  — r 


J>x)  N  (x)  -  J  (x)  N  (mx) 


where  m  *  a/b,  x  *  h2b,  jy  -  h3b  with  y  *|h3b|  .  Using  the  recurrence  relation 
of  the  cylinder  function 


ZnU) 


Z  i ( z)  +  (n/z)  Zn  (z) 


n+1 ' 7  '  “n 

and  the  relation  between  the  Hankel  function  and  the  modified  Bessel  function  of 
the  second  kind 

Kn(z)  ■  eJ^n+1)*  H^(z  ej|,T ) 

eq.  (6-5)  can  be  transformed  Into 
L(y)  ■  R(m,x) 


(6-5) 


where 


K.(y) 

L<y>  ■  y 

and 

J0(x)  N^mx)  -  J^(mx)  NQ(x) 
R(m*x)  =  x  J^x)  ^(mx)  -  J]  (mx)  yx) 


(6-5a) 


(6-5b) 


where  L(y),  the  left-hand  side  function  is  only  a  function  of  y  alone,  whereas 
the  right-hand  side  function,  R(m,x),  is  a  function  of  x  with  m  as  a  parameter, 
x  and  y  are  further  connected  by  the  equation  deduced  from  eqs.  (6- 3a)  and 


(6-3b),  that  Is 

/  •  h2  -  Vo  *  Uy>2  '  Vo 

(6-6) 

hence 

h2  -  (Jy)2  -  (  «r2  -  er3)  k | 

(6-6a) 

Multiplying  both  sides  of  eq.  (6-6a)  by  b2  yields 

x2  +  y2  *  R2 

(6-7) 

where 

fr3 

(6-7a) 

which  together  with  eqs.  (6-5),  completely  determines  the  frequency  and 
parameters  of  the  surface  waveguide.  Notice  that  eq.  (6-7)  Is  a  circle  of 
radius  R  centered  at  the  origin. 

Graphs  of  the  functions,  L(y)  and  R(m,x),  for  four  values  of  m  are  plotted 
In  fig.  (35),  (36),  and  (37),  and  (38).  The  relation  between  the  parameters, 
x  (-  hgb)  and  y  (  “I h J  b),  for  four  values  of  m  are  plotted  in  figs.  (39), 
(40),  (41),  and  (42).  A  graphical  solution  of  the  mode  equations,  (6-5) 
and  eqs.  (6-7)  Is  shown  In  fig.  (43)  for  the  m  -  0.8  case,  in  which  R  has 
arbitrarily  chosen  to  be  26,  this  corresponds  to  b/  ^  =  3.3  for  *  2.56 
and  ■  1.  No  attempt  has  been  made  to  analyze  and  discuss  In  detail  this 
structure,  since  it  has  several  properties  which  are  very  similar  to  the  two 

other  waveguide  structures  discussed  In  sections  IV  and  V.  However,  It  should 
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be  pointed  out  that,  by  inspection  of  fig.  (43),  when  R  lies  between  8.3  <R  ^23.6 
corresponding  to  1.0  £  b/^0  ^  3.0,  no  circles  drawn  with  these  values  of  R 
can  Intersect  any  mode  characteristic  curves  for  the  case  of  m  =  0.8,  and  when 
23.7^  R  <  40  (not  shown  in  fig.  (43),  the  circle  intersects  only  with  the 
second  branch  of  the  mode  characteristic  curve.  This  implies  that  a  single 
mode  operation  is  possible  at  still  much  higher  frequencies  despite  its  short¬ 
comings  discussed  earlier. 

VI.  Conclusion 

It  is  the  objective  of  this  research  study  to  investigate  theorectically 
the  various  ways  and  means  by  which  the  fields  and  waves  generated  from  the 
emitted  electromagnetic  radiation  from  the  Prasad-Leibys'  plasma  tube  can  be 
effectively  confined,  controlled,  and  re-directed  out  the  interaction  volume 
of  the  plasma  tube.  In  view  of  the  fundamental  limitation  of  the  requirement 
of  an  adequate  volume  to  house  the  plasma  tube  no  metal  waveguide  has  been 
considered  due  to  the  co-existence  of  many  modes  in  an  over-sized  metal 
waveguide.  The  waveguide  structures  investigated  belong  to  a  class  of  wave¬ 
guides  which  employ  dielectric  tube  as  their  basic  element,  known  as  the 
dielectric  tube  waveguides.  Three  of  such  waveguides  structures  have  been 
analyzed.  They  are 

Type  I.  The  hollow  dielectric  tube  waveguides 

Type  II.  The  cylindrical  metal  waveguide  partially  filled  with  a 
coaxial  dielectric  tube 

Type  III.  The  hollow  dielectric  waveguide  supported  by  a  conducting 
cylinder. 

In  essence,  Type  II  structure  may  be  considered  as  a  special  case  of  the 
type  I  structure  in  which  the  region  III  is  being  Isolated  by  a  grounded 
conducting  cylindrical  tube,  whereas,  type  III  structure  may  be  considered 
as  a  special  case  of  type  I  structure  in  which  region  I  is  Isolated  from 


regions  II  and  III  by  a  conducting  cylinder. 

The  field  problem  of  the  type  I  structure  has  been  formulated  and  solved, 
and  from  which  the  transcendental  mode  equations  of  the  other  two  structures 
can  be  derived.  These  transcendental  mode  equations  are  solved  graphically 
in  conjunction  with  their  respective  dispersion  equations.  It  is  found  that 
both  types  I  and  II  structures  are  readily  suitable  to  integrate  with  the 
Prasad-Leibys'  tube  to  form  a  complete  millimeter  wave  power  source  unit  much 
like  the  klystron,  or  traveling  wave  tube.  However,  type  III  structure,  in 
view  of  its  structural  limitation,  is  not  suitable  for  this  purpose  unless 
some  structure  alteration  is  made. 

As  a  microwave  or  millimeter  wave  transmission  line  all  three  structures 
are  quite  capable  of  single  mode  operation  at  very  short  wavelengths  when  thi 
dielectric  tube  is  used.  It  is  a  surprise  to  find  that  as  the  operating 
frequency  is  raised  the  metal  waveguide  which  is  partially  filled  with  a  di¬ 
electric  tube,  behaves  like  a  dielectric  tube  waveguide  which  is  limited  on 
the  outside  by  a  conducting  cylindrical  metal  tube,  and  on  the  inside  by  a 
medium  of  smaller  dielectric  constant.  The  fields  and  waves  are  confined  to 
the  space  occupied  by  the  high  dielectric  constant  dielectric  tube. 

Design  information  regarding  to  the  choice  of  the  waveguide  dimensions 
and  the  dielectric  constant  for  the  dielectric  tube  in  relation  with  the 
guide  wavelength  Xg  and  the  wavelength  in  free  space  is  provided  for. 
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